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Matrices having some specified structure are called 


structured matrices, (These are encountered in many fields 
like electrical engineering, structural engineering etc. 

This thesis deals with generation, recognition and storage 
of structured Boolean Matrices, The generation of these uses 
shift -register technique, algorithms for recognising different 
structured matrices from their row or column permuted matrices 
are developed, nlso wo have presented algorithms for recog- 
nizing three of these structured matrices from their row asncL 
column permuted matrices. Another algorithm for reducing the 
Bandwidth of a given random matrix by possible removal of 
rows and columns is included. 

The characterization of classes of genera table 
matrices based purely on their structural aspect is given, 

We -have also given some matrix representations, e»g., 

Polynomial representation, V/alsh-f unction representation otc. 
This was done for the purpose of identifying to which class 
does a matrix belong. 
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CHAPTER 1 


IHPROKJGTIO N 

As is clearly pointed out by the title, the thesis 
deals with goner at ion, recognition and storage of structured 
Boolean matrices* The objective of the thesis is three-fold, 

i) In many application areas, mentioned later, we 
encounter large, structured as well as sparse, matrices. 
Sometimes the storage requirement of these matrices is 
prohibitively large, resulting in our inability to solve 
the given problem, be it important or interesting. 

In such cases, instead of storing the full matrix in 
the core, we cake recourse- to some storage technique, 
thereby reducing the storage requirement of the matrix. Tho 
onco abandoned problem, can now be tackled easily, 

ii) If tho matrix in question is structured, it is highly 
inefficient as far as the execution time is concerned, to 
deal with all elements of tho matrix while performing arith- 
metic operations on it. 

In suoh cases again, we try to reduce bhe computer 
execution time by storing tho matrix utilizing its structural 
property, thereby avoiding as far as possible computations 
with 'zero' elements of the matrix. 

iii) Cfton tho structure of the matrix is not evident at 
the context. In such cases, 


recognition of one of the 
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structures, in tlie matrix is of prime importance. We will 
call this the 'Recognition' problem. 

After o. matrix is recognized to be of one of the 
structured forms, c e wish to recover the structured form 
of the matrix from the seemingly random form. To develop 
methods for doing this was one of our prime objectives. 

Literature survey revealed that though storing and 
indexing techniques exist for sparse matrices, nobody has 
paid any attention to generating, storing or recognizing 
methods for structured matrices. While the number of refer- 
ences involving sparse matrix algorithms are plenty Q 1 , 2, 
3, 8, 9, 11, 14, 18, 30, 52, 35, 41, 42, 43 J ; those 
dealing with structured matrices are only a handful [_ 5 , 15, 
25, 26, 34, 35 j. 

We have considered six structured forms of matrices 

viz. , 

i) Band Diagonal Matrix 

ii) Band Triangular Matrix 

iii) Bloch Diagonal Matrix 
i v ) ' Ov er 1 app ing- blocks' Mat r ix 
v) Mult ip 1 e and 'Overlapping-blocks’ Matrix, and 
vi) Y-shaped Matrix. 

If a matrix is structured one, it is not necessary 
to store the full matrix in the core. In fact, we need only 
one row in the core at a time, and other rows can be 
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generat ed. 

The chapter following the present one gives the 

generating techniques of ’structured matrices’ using shift- 

registers. In this, wo generate (i+1)““ row of the matrix 

from the i — row, using the structural properties of the 

matrix. Then, at any time, storage required is proportional 

to N, while if the full mo.tr ix were stored; storage required 

2 

is proportioned, to N . Algorithms for generating all the 
above mentioned structured matrices arc developed and 
presented. Also, it is possible to recover all of them from 
their row or column permuted matrix. Algorithms for doing 
this, are also given. 

It has boon observed that sometimes both row and 
column permutations are necessary for reducing a given 
matrix to one of tho structured forms. Chop l or 3 gives 
three such algorithms to obtain, 

i) Band Diagonal 

ii) Band Tri . ngular , and 

iii ) ’Ovcrlapping-Blocks ‘ , 

matrices from their row and column permuted matrix. 

In some cases, it is not possible to reduce the 

matrix into either of iihe structured forms by row and 
column permutations. In suoh oases, a few rows and/or 
columns are to bo removed and the remaining matrix converted 
to one of the structured forms. An algorithm for such 
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conversion to Block diagonal matrix exists and 

a new algorithm for conversion to Band diagonal matrix is 
given in Chapter 4 « 

It is poo si lie to generate some classes of matrices 
using the generating scheme of shift registers and a few 
allowed operations. The description of these classes, based 
purely on their structural aspects, appears in Chapter 5* 

To charac Leri 550 these classes mathematically? we 
tried some representations of matrices, e»g», Polynomial 
r ©present at io n, Walsh-function representation etc. ‘Those 
representations are discussed in Chapter 6. It should how- 
ever be mentioned that these proved to be fruitless olid a 
method of characterizing the classes of matrices considered 
in this thesis, remains unsolved. 

The complete work is summarized in Chapter 7? which 
is the concluding chap tier. 

To quote Tewarson j]5l] , sparse matrices occur in the 
solution of many important practical problems, o.g., in struc- 
tural analyses, network theory and power distribution systems, 
numerical solution of differential equations, graph theory, 
reactor diffusion as well as in genetic theory, behavioral 
and social sciences, and in computer programming. As the 
technology increases in complexity, wo can expect that large 
sparse matrices' will continue to occur in many future appli- 
cations ' involving large systems, e,g», simulation of traffic 
lights, pattern recognition, urban planning and so on. 
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According to Pooch and Jtfieder , {* 2lJ , it has been 

observed that in nany of the application or eas mentioned 
above, the matrices are not exactly random ones. In elec- 
trical network and power distribution systems, the matrix 
is generally in random bond or block- diagonal form with 
the elements representing circuit voltages, impedances or 
power sources [pS] j in structural engineering, the sparse 
matrix is generally of band or block form, with the band- 
width or block dimension representing the number of joints 
per floor [j 3, 19] » in reactor diffusion problems and diffe- 
rential equations, the band form of matrix is most common, 
with the bandwidth being the number of points used in a 
point-diff ercnc o formula [23, 24,37j. 

o, g. , If wc solve the boundary-value diff orential equation 
2 + 3t 2 = y+yVy" by using a central difference approximation 
(3 points), and 10 intervals between the points y(t-O) = 0, 
and y (t — 1 ) = 1; tho resulting 9x9 matrix obtained is shown 
on page 6, It is a Band Diagonal matrix. 

A 5- point interpolation would yield a bandwidth of 
5,' 50 int envois would result in a 49^49 matrix. 

In favour of Band Diagonal matrices, it can be said 
that in most oases, they are processed either wholly by 
rows or columns, -and little or no processing of single 
element occurs. Hence at any time, only a few rows/ columns 
need be maintained in fast core for immediate use. 
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-199 100,5 


99 -199 101 

* 

98.5 -199 101.5 

98 -199 102 

97.5 -199 102.5 

97 -199 103 

C \ ' 96.5 -199 

96 

I 

» 

I 


103.5 
199 104 

95.5 -199 


Thcrofoxv., it is advisable that whenever a user 
encounters such a matrix, a special effort is desired on 
his part to add IP t his processing algorithms to the case at 
hand. This adaptation should bo made because of tho inherent 
simplicity of processing, manipulating and solving band 
matrices, and also because of the opportunity to minimize 
core alio cat ion and execution time. 

a,s already seen, most of the structured matrices 
that one encounters in different application areas arc 
sparse, too, This is especially true in the case of large 
matrices because of the fact that the matrix density 
(defined as the number of nonzero elements of the matrix 
divided by the total number of elements in the full matrix) 
decreases as matrix size increases. 



7 


Now, if number of ’zero' entries in the structured 
matrix, within a band or a block as the case may bo, is 
quite large, then it is not advantageous to store the matrix 
using its structural property. It pays to treat such matrices 
as sparse matrices only, and use the existing indexing 
techniques for storing sparse- matrices. 

To moke this report self-explanatory and self- 
sufficient to the reader, a brief account of existing index- 
ing techniques to store sparse matrices is included hero, 

The different techniques are explained below [21], 

1 • Bit -Map Sch eme: - 

£ 

The given matrix is first converted to a Boolean 

matrix by putriny nonzero entries as 1’s. These elements 
arc stored in an array. A row-index vector is created to give 
address of the first nonzero entry of that row. 

e.g. 
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R ow- Index A ctual 

Vector Values 


The rows of hit nap arc packed to store maximum, 
number of rows in a word. 


e.g. , 


10 10 0 10 1 


1 0 0 1 0 0 1 0 


Bytc^ Byte Byte 

1 2 3 


If I and J are dimensions of the matrix, end B is 
the number of bits ./word, then this scheme requires (I*J)/B 
= W words. Storage required for row-index vector = I*A/B 
words, whore A - number of bits required for an address. 


. *, Percentage storage required of the full matrix 


100 

~B 



Using B = 32 bits and A = 16 bits for a. 500 x 500 
matrix, 5*/* dense,* total storage required is 12,500 (for 
nonzero elements) + 8313 (for bit-map and row-index vector) 
(= 20813 words - 8.325*/* of the full matrix. 
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If we wish, to access a- • , we note that this hit 

J. J 

lies in S. =5 (i-l) * J + j + (B-1) 2 /B word, of the hit 
nap. The required bit can now be isolated by masking all 
other bits by a logical operation or shifting the word. 

If we store the bit -nap by rows, then to perform 
a column operation, we shall have to count j — bit for all 
I rows, which is a time consuming process, especially if 
such operations arc quite frequent. 

A dvantages 

1. I'his scheme of indexing requires a definite core 
storage. 

2. Row-access tine off erred is very reasonable. 

3. Quito core efficient if matrix density is greater 
than 5/. . 

4. If operations to be performed are mostly row. 
operations, then it is quite fast as regards execution time. 

5. fast Boolean operations arc possible. 

U If stored by rows, it has a very poor column 
access time, 

2. Reordering data elements requires considerable time. 

3. Makes poor use of parallel processing, 

4. If matrix density is very low, the scheme does 
not prove to be core efficient. 
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II* , Address-Map Scheme : - 

This indexing scheme is quite similar to one already 
described; the difference lies in the fact that here, we 
store an address or an address displacement for each matrix 
element. Obviously, this requires 0 times more core to 
store the matrix, whore 0 - number of bits used for an 
address or address displacement. 

If there are less than 64 nonzero entries per row, 
the address displacement requires only 6 bits. 

The overall percentage storage requirement of the 
full matrix in this scheme with the row-index vector is 


q — J.22 r o + */ 

add- map B L ^ + j J A • 
th 


To retrieve a. clement, it is necessary to access 
the s i 5 (A- 1 ) * J + d j 


character or byte (assuming 

one byte used for an address displacement). If the entry 
"bid 

in S. -- byte is sort', the matrix element was null; other- 

"t ll 

wise the constant of — byte is added to the row-index 
element to give the address of the nonzero element. 


JX dv aht ages ;- 

1. In machines where character addressing is available, 
the scheme mahos very efficient use of parallel processors. 

2. Any element of the matrix can be readily accessed, 

3. As compered to hit-map scheme, it requires less 
execution time to reorder data elements. 
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4, Row and column access times exhibited are very 
reasonable. 

D isadvantages 

1« It requires a large fast core for indexing, 

2, As compared to 'threaded-last* scheme, it requires 
large execution time to reorder matrix elements* 

II I , Delta or Displacement Indexin g 

This indexing technique is almost identical to 
Address-map scheme. Instead of storing absolute addresses 
of matrix elements, address displacements from previous 
one are stored. In one of the existing delta-indexing 
schemes, one b.' r -bit extended precision word contains one 
16-bit index and six 8-bit displacements to the index. 
Hence column indices of 7 elements can be referred to by 
loading and processing one extended precision word. 


e.g. 

773 797 806 873 

col, index 1 0,1, 2 0.1,3 0.1,4 

Blocked Index Wor d 

The corresponding delta-indexed word is 
773 24 9 67 

0*1,1 delta delta delta delta 


The two indexed words are equivalent. 
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I Y . T hread eel -List Scheme ; - 

This indexing method concentrates only on the nonzero 
elements of the matrix. Lor each of these, an element of 
an array in core is provided. Each such element consists 
of at least 3 components; one containing row and column 
indices of the element, second containing its value and the 
third pointing to the next array element, 

e.g. ? 


Address 

row 

col. 

value 

po inter 

773 


31 

- 

00 

o 

VO 

I 

t 

- - - - 

- - - ~ 

- - - - 

t 

- J 

1 

t— -► 809 

27 

3 

- 

976 




i 

t 

I p* 883 

33 

25 

JZ 1 

- 

983 

1 ; 

T t 




f 

1 r_ _ _ _ _ 


- - - - 



T 

U- 976 

4-1 

19 


883 


r 

t 

t. 


983 






Minimum storage required fur indexing in threaded- 
list scheme is 

Y/ = £ (J * T + A) D + Y 3 words; 

where I, J are dimensions of the matrix 

A = Number of bits required for an address, 

B t= Number of bits/word, 

T 5 = Number of bits used for column index element, 

V - Number of bits used for row index clement , 

D *= Density of the matrix. 

Advantages 

1, This is the sole technique that allows a simple 
and fast execution method of reordering; adding or annihi- 
1 at i ng dat a eL em ent s , 

Pi s advant a/ ; eg : - 

1 . It requires large core, especially if the matrix 
density is high • 

2, Slower access bine for columns as compared with 
that in r row- column 1 scheme, 

Y , Row -Column Schem e ; - 

There exist plenty of row-column schemes differing 
from one another in construction of different arrays {jf]. 
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Scheme 1: 


To each nonzero element of the matrix, there corres- 
ponds an item of two storage cells 5 the first containing 
the row-index, and the second, the value of the element. If 
first element is r 0 f , we interpret it as end of the current 
column, in which caso the second element indicates the index 
of the next column. Zeros in both cells indicates end of the 
matrix storage. 


entries 


e 


rr 


• fa* 




If n = dimension of the matrix, p = no, of nonzero 
then storage required «= 2 (n + p + 1) cells. 

PS 0 3 0 | 


A — 


0 

4 


2 

0 

0 


0 

0 

8 


7 

6 

0 


For *1, storage goes as follows. 

(0,1; 1,5; 3,4; 0 , 2 ; 2 , 2 ; 0,3; 1,3; 4 , 8 ; 


0,4; 2,7; 3,6; 0 , 0 ). 


Scheme 2 : 

Three arrays, namely, YE (Value of elements), 

RI (Row indices) and CIP (Column index pointer) are created. 
Rl(a) contains row-index of the corresponding element YE (a) 
of YE. OIP indicates position of first element of a column 
in YE array. 


Total storage required is clearly 2p + n cells 
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For A, 

• VE *= (5, 4, 2, 3, 8, 7, 6) 

RI = (1, 3, 2, 1, 4? 2, 3) 

CIP = (1, 3, 4, 6) 

S ohene 3 : 

In this scheme, with each nonzero element of the 
given matrix, a unique integer 0. . is associated. 

XJ 

0 ij = i + (o " 1) n. 

Two arrays VE ...nd TH (Th of a) are created and stored. 

for a, 

VE = (5, 4, 2, 3, 8, 7, 6) 

TH = (1, 3, 6, 9, 12, 14, 15) 

To recover the original matrix; j is the least 
integer > 0(i,j)/n and i = 0(i,j) - (j-l) n. 

To store a, matrix in this' case, we require only 2(3 

cells . 

Scheme 4 : 

This indexing soli an e can ho used if the matrix is 
symmetrical. In that case, only the lower triangular pant 
of the matrix is stored. The storage consists of two .arrays, 
VE and PD (Position of Diagonal elements in VE) . Por ca,ch 
row, the leftmost nonzero element and all other elements to 
its right upto and including the diagonal element are stored 
in VE. 
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Thus the i~~ row of the matrix needs (9-t-l) locations, 
n 

and YE will have S (9,-+ 1 ) elements j where 0.- = distance of 

i=1 

x th 

leftmost nonzero element from its diagonal, in i — row* 


n 

To store the matrix then, requires 2 9 . + 2n loca- 


tions. 


e.g., 


Then, 


(4, 7, 3, 0, 5, 2, 0, 1) 
0 , 2, 5, 8) 


An element a, . can be recovered as follows. 

■p-L J 

PD(i) - (i-j) is the position of a. . in YE jjrovided 

-L J 

PD(i) - {i-j ) > PE (i-1); i, c. does not lie to the left 
of the first nonzero element of the i — row, otherwise* 

a iD = °* 

e.g., To recover a ^ from YE, we have, 

PD(4) - (4-3) « 8-1 - 7 > 5 = PD(3) 

, , is in YE (7 ) . 

The discussion of different indexing techniques indica- 
tes that an exact comparison of execution times must reflect 
the type of mathematical manipulation being performed on the 



17 


sparse matrix, c . p . , The- bit-map method is of particular 
use when the matrix is used to produce an ’’optimal" ordering, 
so the matrix inverse will not have a great increased density. 
In contrast , the row-column method is faster than other 
methods when manipulations involve onu row (col.) acting on 
other rows (cols.). 

The second important aspect of indexing scheme 
Selection is the conservation of execution time. If arith- 
metic operations are to be performed on the data, primary 
consideration should t ,> given to a row-column method; if 
Boolean arithmetic or reordering algorithms are to bo performed, 
the bit map scheme should be considered first; and if a great 
number of data elements are to be reordered, created or 
annihilated, a threaded-list scheme deserves first considera- 


tion. 


CHAPTER 2 


GENERATION OE STRUCTURED MATRICES 


Given a random matrix, wo first convert it to a 
Boolean matrix by putting 1 in place of each nonzero entry 
of the matrix. If our generation procedure is row-wise, we 
store an array VE (Value of Elements) which contains actual 
value of nonzero elements appearing in row- wise order. 

fch 

When we wish to perform operation on i — row elements 

"ij i t 

we successively generate first, second, (i-1) — row 

keeping count of total number of nonzero elements occurring 

"fc ti 

in these rows, let K be this number. Then we generate i — 
row and obtain actual values of nonzero entries from VE(K+1 ) 
onwards . 

Thus, extracting uxact values of matrix elements from 

Boolean rnatiix representation can be done s imply by keeping 

count of total nonzero entries in previous rows and then 

\ 

retrieving elements from th > array, VE, Keeping this in mind 
henceforth, we shall be confining our attention to Boolean 
matrices only. 

Figure 1 shows a few "s ci ueturu', -matrices' 1 . 



^- s ) - JB1 ocA i 
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It can be argued that not all of these occur in 
applic at ion- areas . However, as we shall consider later 
(Chapter 6), the combination of two of these matrices with 
AID, OR or EXCLUSIVE OR element -wise operation between them 
yield seemingly random matrix . This leads us to storing two 
structured matrices instead of a random matrix, which in most 
of the cases involving large matrices , results in reducing 
core-storage requirement . 

The generation of those structured matrices uses Shift- 

S "t 

register technique; (i*1) — row being generated from the i — 
row using the structural property of the matrix. The generating 
scheme is shown in Fig #2, 



FIGURE 2 
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The operations allowed are as follows :- 

1 ) Shift right 

2) Shift left 

3 ) Count and 

4 ) Compare . 

The importance of working with structured matrices 
is quite evident. As these can be generated row-by-row by 
using shiit registers, it is not necessary to store the 
complete matrix in core at any time. Only the pertinent row 
can be generated and used for further computations. Clearly, 
this means an increase in Execution time, only as fsr as gene- 
ration of consecutive rows is concerned. However, a problem 
which could nob be solved earlier due to prohibitive demands 
on memory requirements, can now be taken up quite conveniently. 

If a given sparse Boolean matrix is one of the struc- 
tured-matric.es, or can be converted to it by row-column permu- 
tations, (algorithms for a few random matrices are given in 
this and the next chapters); then the total storage required is 
proportional to M , the size of the matrix, while if full 

2 

matrix is stored, the storage required is of the order of N , 
which is considerably large, especially in cases when N is 
quite large. 
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(i) 

Step 

Step 


Step 

Step 

Step 


ALGORITHMS FOR GENERATION OP STRUCTURED 
MATRICES USING SHIFT REGISTER TECHNIQUES 

I 

Band-Diagonal Matrix : 

Is Given I,J,N. Initialize a count K=1. 

2; Set the shift register contents to I 1 *s 
and (N-l) 0’s; i.e., the first row. 

3: PRIED the shift register contents. 

4: If all rows over, GO TO step 9- 

5; Increment K by 1 and shift right the contents of 
the shift register. 



Step 6s If Kf>J, GO TO Step 8. 

Step 7i Append a " 1 n in first position of shift register. 


Step 8; 


GO TO Step 3. 

Append a "0" in first position of shift register. 


GO TO Step 3 * 
Step 9: STOP. 

e .g. , , I = 2 , J * 3 , 

Pirst row is 
The next two rows are 
Then, the remaining 
rows have "0" in first 
bit position 


I = 6 

1 

1 

1 

0 

0 

0 


1 

1 

1 

1 

0 

0 


0 

1 

1 

1 

1 

0 


0 0 

0 0 
1 0 
1 1 

1 1 

1 1 


o" 

0 

0 

0 

1 

1 


*Read "Band Diagonal Matrix" as "fully-packed Band Diagonal : 

.Matrix” if Bands are all I's. j 

i 

i ! 


1 
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The matrix generated is the desired Band-Diagonal 


matrix . 

(11) Block-Disgunal Matrix 


*1 


Step Is 

Step 2% 

Step 3: 
Step At 

Step 5: 


Step 6: 
Step 7: 


Step 8: 
Stop 9: 

e .g . > 


Given the size of the matrix N 
and block sizes. 


/ 


Ji : 


d J ^ 










Set the shift register c< ntents to first row of the 
first block. PRIM) it. 

If all rows over, GO TO Step 9. 

I± the present block generation is coinplecc, GO 
TO Step 6. 

PRINT the shift register contents. Keep the count 
of number of rows generated of the present block. 

GO TO Step 3. 

If all "blocks 11 generated, GO TO Step 8. 

Set the shift register contents t > first row of the 
next block. PRIM? the shift register contents and 
GO TO Step 3* 

Generate the remaining "2 jER 0" rows , if required. 
STOP. 


i-j = 2. 


U = 3. 


i 2 =1. 


= 2 ' 


i, = 2. 
3 

i™ = 1. 


N = 7. 


'"1 
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it row of the first 
hk is 

5b such rows form first 
he 

Nbluck is generated 

O Oil* 


(■ Band-Triangular Matrix : 
St; Given triangle sizes 


1 10*000 

1 1 0 0 0 0 0 

1 1 0 0 0 0 0 

0 0 1 0 0 0 0 

0 0 1 0 0 0 0 

0 0 0 1 1 o o 

0 0 0 0 0 0 0 

J [ ~ 0 




dimension of the matrix. 


St; 



first hit position and all 0's, u ., the first row 
of the first triangle. 


St; PhlNT the shift register content 
St; If all rows over, GO 10 Step 10. 

St; If all rows of the present triaie over, GO TO 
, Step 7. 

Sts Shift right the shift register Cents , and put 
a "1” in the column in which thowas first "1" 
in the last shift register corite. GO TO Step 
St.e If all triangles generated, GO Ttep 9. 

Sts Set the shift register contents first row of 
• the next triangle, and GO TO Ste. 


f 


l 
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Step 9: Generate the remaining "ZERO' 1 rows, if required . 

Step 10; STOP. 

e.g., J(1) = 0 J(2) = 3 J(3) = 2 J(4) =2 N = 7 

First row is 1000000 

Then the first 1 1 0 0 0 0 0 

triangle is 

generated 1110000 

Then the next 0001000 

triangle is gene- 0001 100 

rated and su on, 0 0^ 0 0 0 1 O 

0 0 0 0 0 1 1 



Step 1: Given number of sets of blocks, starting and 

finishing column of each block, and size of the 
matrix . 

Step 2* Set the shift register contents to first row using 
"starting and finishing columns of blocks" 
information. 
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Step 3: PRINT the shilt register contents. 

Step 4: If all rows over, GO TO Step 10. 

Step 5: If Che present set of blocks over, GO TO Step 7. 
Step 6; After keeping count of tho number of rows of the 
present set of blocks generated, GO TO Step 3. 

Step 7; If all sets uf blocks over, GO TO Step 9« 

Step 6: Set the shift register to first row next set of 
blocks, age in by using "storting and finishing 
columns of each block" inf .mnation. GO TO Step 3« 
Step 9: Generate the remaining "ZERO" rows, if required. 
Step 10: STOP. 


e.g., ^11 

= 1 

f 11 

3 


S 12 “ 

5 


M3 

11 

CM 

<H 

S 21 

= 2 

f 21 

4 





S 3 1 

= 2 

II 

K\ 

4 


s 32 = 

6 

f 32 ~ 7 

h : 

= ^2 = 

s = 2 

• 


II 

• 


Then, thi 

e matrix 

generated 

is 





1 

1 

1 0 


1 

i 


0 



i 

1 

1 

1 0 


1 

i 


0 



0 

1 

1 1 


0 

0 


0 



0 

1 

1 1 


0 

0 


0 



0 

1 

1 1 


0 

1 


1 



0 

1 

1 1 


0 

1 


1 


1 

0 

0 

0 0 


0 

0 


0 
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(v) ’’ Overlapping-blocks 11 Matrix : Same as (iv) 

(vi) V-Shaped. Matrix: 



Step 

1: 

Given I, J, R. 

Step 

2: 

Set SRI to I 1 r 



0 f s and I 1 1 s . 

Step 

3: 

Concatenate SRI 

Step 

4: 

If all rows ove 

Step 

5! 

Increment K by 


that of SR2 left, dropping the rightmost bit of SRI 
and leftmost bit of SR2. 

Step 6; If K>J, GO 10 Stop 8. 

Step 7: Append a ”1” in first bit position of SRI and last 
bit position of SR2. GO TO Step 3. 

Step 8; Append a "0" in first bit position of SRI and last 
bit position of SR2. GO TO Step 3. 

Step 9: STOP. 


e.g. , 


With 1=1, J = 3 and N = 8; we get the matrix 


1 
1 

1 1 
0 1 
0 0 
0 0 0 


0 0 0 0 
10 0 0 
1 0 0 
1 1 1 
1 1 1 
1 


0 0 


0 1 
1 1 
1 1 
1 0 0 
10 0 0 


1 

1 

1 

0 


00000000 

[oooooooo^ 
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If the given matrix is net in one of the structured 
forms but is row (or column) permuted matrix of one of such 
forms, then it is always possible to recover the structures 
matrix from the given matrix by row (jr column) transformations. 
The algorithm fur doing the same are presented below. 

RECOVERY OH STRUCT DRIP MATRICES BY ORLY 
ROW PERMUTATIONS OR A GIVEN MATRIX 

( i ) Band Diagonal Matrix : 

Step 1: Given M. If the column codes of all rows arc not 
continuous GO TO Step 21. Otherwise, calculate 
the row-sums and let MAXR = max. of the&e. 

Step 2; Consider the rows having 1 in the first column. 

Of these row-sums, let K bo the minimum. 

Step 3: Rind a row whose column codes are 
If not found, GO TO Step 21. 

Step 4: Put this row in sorteu row (SR) vector. If R^MAXR, 

GO TO Step 6. 

Step Si Increment K by 1 and GO TO Step 3« 

Step 6: If MAXR = I, GO TO Step 16. 

Step 7s 1 = 2; J = 1. 

Stop 8: Rind a row having column codes I, 1+1, . . . ,(MAXR+J) . 

If not found, GO TO Step 21. 

Step 9: Put the row in SR. Increment J by 1. If all rows 
over GO TO Step 20. 

Step 10: If (MAXR+J)>N, GO TO Step 12. 
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Step ,1: Increment I by 1 anu GO TO Step 8. 

Step 12: Assign (N-MAXR+1) to K. 

Step 13: Find a row having column codes K, K+1 , 5T. 

If not fuund , GO TO Stop 21. 

Step 14: Put this row in SR. If all rows over, GO TO Step 20. 

Step 15: Increment K by 1 and GO TO Step 13. 

Step 16: If there exists a row having column codes 1,2,...N; 
GO TO Step 18. 

Step 17: II = 2. GO TO Step 13. 

Step 18: Put the row in SR. If all rows over, GO TO 20. 

Step 19: GO TO Step 16. 

Step 20: PRINT the row- permuted matrix and STOP. 

Step 21: Report that given matrix can' r t be converted t_> 
Band-Diagonal form, and STOP. 

( ii ) Block Diagonal Matrix : 

Step 1: Given M. If column codes of all raws aro nob conti- 

nuous , GO TO Step 7. Otherwise, calculate raw-sums 
and put rows having ZbRO row-sum ac end posiriuns 
of the sorted row (SR) vector. 1=1. 

Step 2: Find a raw having column codes as I, 1+1,.... If 

no such raw found, GO TO Step 7» 

Step 3: Put this row IR in SR. Ii all rows over, GO TO 

Step 6. 

Step 4: Find a row having zero distance (i.e. identical row) 

from IR. If such a row exists, GO TO Step 3. 
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Step Let highest column, cole of IR be J. Put I = J+1 

and GO TO Step 2. 

Step 6: Print the permuted matrix and STOP. 

Step 7: Report that matrix can't be converted to GuSire.d 

form and STOP. 

(iii) Band Triangular Matrix : 

Step 1; Given M. Place "ZLRO" rows, if any, at the end 
positions of the sorted row vector (SR). 

1=1 J = 1 . 

Step 2: Find a row having column codes I, 1+1, . . . , J. If 

no such row found, GO TO Step 5. 

Step 3: Put this row in SR. If all rows over, GO TO Step 8. 

Step 4: Increment J by 1 and GO TO Step 2. 

Step 5: If (J-I) = 0, GO TO Step 7. 

Step 6; I = J and GO TO Step 2. 

Stop 7: Report that matrix is at convertible to uooired 

form and STOP. 

Step 8: PRIM' the ruw permute J matrix and STOP. 

( iv) 21 

Step Is 



Step 2: 


I ,= 1 . 
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Step 3 } 


Step 4: 

Step 5: 

Step 6: 


Step 7: 


Step 8; 

Step 9: 


Find a row having column codes I, 1+1,.,. . In case 
of clash, choose the row having minimum row-sum. 

If no such row found, GO TO Step 8. 

Put this row in Sorted How Vector (SR). Call it 
IR. If all rows over, GO TO Step 9* 

Pincl a row having zeru distance from IR. If such a 
row exists, GO TO Step 4. 

Pina a row which shares one or more Column with IR. 

In case of clash, seluct one having maximum columns 
in common with IR. If clashes occur again, choose row 
with minimum row-sum. If such a row exists, GO TO 
Step 4. 

Let J be the highest column code of IR . I = J + 1 
arid GO TO Step 3 . 

Report that matrix cannot be converted to desire>d 
form and STOP. 

PRINT the row permuted matrix and STOP. 


(v) " Multiple and Overlapping -blocks Matrix " j 

Similar t j ( iv) . 


(vi) V-Shaped Matrix : 

Step 1: Given M. Calculate r^w-sums . Place "ZERO" rows in 

the end positions of sorted row vector (SR). N = size 
of the matrix. 1=1. 
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Step 2: 

Step 3: 

Step 4} 
Step 5: 
Step 6; 

Step 7: 
Step 8; 
Step 9: 
Step 10*. 
Step 11: 

Step 12: 

Step 13: 
Step 14: 

Step 15: 


Find a row having column codes I, 1+1, ...I+K, 

N-K, . 11 with minimum K. If nj such row exists, 
GO TO Step 14. 

Induce the row in SR. If all rows over, GO TO 
Step 15* 

Increment K by 1. If I+K <. |~ -j| j ; GO TO Step 6. 

Set K= I"!]- I. 

Find a row having column codes I, I+K, N-K, 

*..,K. If such row r exists, GO TO Step 3- 
Increment I by 1 . Decrement K by 1. Set J = 1. 

If (I+K) <T|1 , GO TO Step 11. 

If I>|" f 1 > S° 10 Step H. 

Se tK=rf-]-i. 

Find row having column codes I, I+K, N-K- J , 

N-J. If no such row exists, GO TO Stop 14. 
Include this row in SR. If all rows over, GO TO 
Step 15* 

Increment I & J by 1 anc GO TO Stop 8. 

Report the impossibility of converting the matrix 
in ecsired form and STOP. 

PRllJT the row-permutd matrix and STOP. 


Tewarson |^3 ij has given a few more desired forms of 
the matrices, which are reproduced, in Fig. 3. 

Given the band, block and border dimensions, it is 


evident that these dorms can be easily generated using shilt 
register technique. 
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CHAPTER 3 

RECOGNITION OF STRUCTURED MATRICES 


It has been seen that the generation, and recovery 
of structured matrices, by row-permutations only, is possible. 
Unfortunately, the matrices that one encounters in real-life 
problems do not always appear m such a concise form. In that 
case, the given matrix is to be Converted, to one of the 
structured forms , b;y employing both, ruw and column permuta- 
tions . The algorithm fur u^ing this are developed fjr three 
structured matrices, viz., 

i) Band Diagonal Matrix 
li) Band Triangular Matrix 
ill) "Overlapping -Blocks" Matrix 

Such recovery of structured! matrices is of great 
importance. Bor example, when solving a large system of 
equations, if it is possible to reduce the coefficient 
matrix into a block diagonal form by row-coiumn permutations, 
the blocks formed clearly inuicatc which equati -ns are to 
be solved for finding which unknown. The technique is very 
well developed and explained by Steward j~27 , 2sQ . Similar 
advantages are noticed for Band Diagonal matrix which appears 
quite frequently in Electrical Engineering and control 
problems, as well as some of the ether application areas. 
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ALGORITHMS FOR RECOVERY OF STRUCTURED MATRICES 
PROM A RANDOM MATRIX, BY ROW-COLUMN PERMUTATIONS 

( x) Bfand Diagonal Matrix : 

All matrices, which can be converted t j Band 
Diagonal form, necessarily satisfy conclitijn C. 

Qjnditiun 0 ; (Does not apply for MAXR or MAXC = N) . 

If MAXR = maximum row sum, then number of rjws 
having row-sum = MAXR shjuld be (B-MaXR+ 1). Similarly, number 
of columns having column sura = MaXC should bo (R-MnXC+l), 
where MAXC = maximum of column sums. The remaining rows/ 
columns should have proper row-s urns/ column- sums so that these 
rows/columns can be placed properly. (i,e., Number )f rows 
having- row-sum = MAXR-1 shwuld be one or two . If ir is two , 
rows having row— sum = MAXR-2 are lound , and Sl. on. If at any 
stage, number uf rows found = 1, then remaining rows should 
have row-sums = MAXR-K , MAXR-(K+1) and so on). 

Step 1 ; Given the matrix M. Compute its row-sums and 
coluinn-sums . 

MAXR cs maximum uf row- sums, 

MIBR = minimum of rovy-sums . 

MnXO = maximum of column sums . 

MING — minimum of column sums. 

Step 2: If condition C not satisiied, GO TO Step 19* 

Step 3s Select a roy v having minimum row-sums. Gall it R1 . 

Put it in sorted row vector (SR). Set IP to MINR. 
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Step 4: Find a row having row-sum = IP + 1, and which is 
nut yet considered. Call it R2. 

Step 5: If R2 is not found, GO TO Step 19* 

Step 6: If distance between Rl and R2 ^ 1 , GO TO Step 4. 

Step 7} Let Rl = R2. Put Rl in SR. Increment IP by 1. If 
IP<;MAXR, GO TO Step 4. 

Step 8; Find a rjw having row-sum = MAXR and distance from 

Rl = 0 (in case when MAXR = N, only) or 2, If no 

such row found, GO TO Step 19* 

Slop 9: Call the row Rl mid put it in SR. If all rows over, 

GO TO Step 14. 

Step 10: If number of rows having row-sum = MAXR already 

considered is not greater than N-MAXR, GO TO Step 8. 

Stop 11: Set IP to MAXIM. 

Step 12: Select a row having r iw-sum = IP, which is nit yet 
considered, and has distance 1 frurn Rl . If n^ such 
low found, GO TO Step 19. 

Step 13: Call it Rl and put it m SR. If all rows aro not 

considered, GO TO Step 12, after decrementing IP by 1. 
; 

Step 14 i 11 columns of M are considered, GO TO Step 1b. 

Step 13*. Transpose M to get M* . IED = MAXR - MINR.Mi.XR = MAXG 
etc* GO TO Step 3- 

Step 16: If MI NO = IRD+1 , GO TO Step 18. 

Step 17: Reverse the column order. 

Step 18: PRINT' the permuted matrix and STOP. 

Step 19: Report that matrix is not c invertible to Band 


Diagonal form and STOP. 
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(ii) Band Triangular Matrix ; 


Step 

Is 

Step 

2: 

Step 

3: 

Step 

4: 

Step 

5: 

Step 

6: 

Step 

7: 

Step 

8: 

Step 

9: 

Step 

10: 


Step 11: 


Step 12; 


Given the matrix M. Calculate its r ;w-sums. 

Put rows having "ZERO" row-sum at the end of 
sorted row vector (SR). Set INDEX = 0. 

Find a row whose row-sum = 1 and which is not yet 
considered. If no such row exists, GO TO Step 12. 

Find the column having non-zero entry. Call it J. 

If the column ,.s already considered} GO TO Step 3* 
Increment INDEX by 1 . Put selected row IR in SR 
vector. Interchange columns J and INDEX. 

If all rows over, GO TO Step'll* 

Set K = 2. 

Find a row having row-sura = K and which is not yet 
considered. If no such row found, GO TO Step 3* 

Call new found row, I. If distance between rows I 
and IR is not equal to 1 , GO TO Step 8. 

Find the column J such that M(I,J) = 1 but M(IR,J)^1. 
Increment INDEX by 1 anu interchange columns J and 
INDEX, Put I— row in SR. Let IR = I and K = K+1. 

If all rows not over, GO TO Step 8. 

Columns are alr&dy permuted by interchanges. Permute 
rows as given by sorted row vector and PRINT the 
permuted matrix. STOP 

Report that matrix is not convertible to Band 
Triangular form and STOP. 
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(in) " Overlapping-Blocks" Matrix ; 


Step Is 

Step 2s 

Step 3s 

Step 4$ 
Step 5i 

Step 6: 

Step 7: 

Step 8: 

Step 9s 


Step 10: 


Given the matrix M. Calculate row-sums and 
column-sums , 

Form blocks of rows and columns and lot BLR(I,J) 
and BLC(I,K) give rows and columns respectively 

i” Vi 

which form I— block. 

Pur eech block, calculate the number of blocks 
which have at least one column in common with it. 
Let KT(I) give this information. 

If for any of the blocks, KT(I)^>2, GO TO Step 10. 
Choose a block whose KT(.) = 0 or 1 i.e. indepen- 
dent block or sharing at least one column with one 
block only. If no such block found, GO TO Step 10. 
Call the block IB and put it in sorted block vector 
(SB). If all blocks over, GO TO Step 9- 
Find a block which has £t least a column m common 
with block IB. If no such block found, GO TO Step 5 
Put the new block IB m sorted block vector. If all 
blocks not over, GO To Step 7. 

SB vector gives blocks in order in which they would 
appear in the final matrix. Permute columns so that 
shared columns bet ween blocks come at proper places 
thus forming sorted column (SC) vector, Form sorted 
row (SB.) vector of sorted rows from BLR. Add "ZERO" 
rows and/or columns if necessary, and print the 
permuted matrix. STOP. 

Report that matrix is not convertible to desired 
form and STOP. 



CHAPTER 4 

ALGORITHM FOR REDUCTION OE L GIVEN MATRIX 
TO BAND DIAGONAL POEM 


In the previous chapter, we presented algorithms 
to convert a random matrix to one of three structure! forms 
by row and column permutations. However, it is obvious that 
we cannot always reeuce a given random matrix to one of the 
structured matrix forms. This can be mainly due to two 
reasons . 

i) Our constraint that all elements of the matrix within a 
block or a band, as the case may bo, are nonzero. 

ii) The matrix itself is random. 

So now, we relax our condition and say that the 
elements in a band or a block can be either 0 or 1, but the 
elements outside it are O's. Even with the constraint relaxed, 
it is quite possible that a matrix cannot be converted to 
desired structured form. In that case, we have to remove a 
few rows and columns, place them as '’borders", and convert 
the remaining matrix to one of the desired forms. 

Rajaraman and Muthukrishnan £ 22j and Weil and 
Kettlerj^’ 39 ,40j have developed algorithms to convert a ranuom 
matrix to block diagonal form by possible removal of a few 
rows which form a border, (figure l). 

A new algorithm is developed to convert a random 
matrix in either a Band Diagonal form (if possible), or a 
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Figure. 1 





singly or doubly bordered band form (figure 2); the latter 
requiring removal of rows or/ani columns to reduce the band- 
width of the matrix. 

fs the nonzero elements of the matrix are now concen- 
trated in a narrow band centered along tho principal diagonal, 
the representation of the matrix on e computer or on paper 
is concise and the arithmetic entailed in the computations 
with it is considerably simplified. 

In the algorithm, we have limited bandwidth to the 
maximum of row-sums and column-sums. If it increases this 
number, rows or/and columns are removed to decrease it. The 
following definitions and criterion will help understand the 
algorithm better. 
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jDefinit ions ; 

k row/column is said to be adjacent to another 
row/column it there exists a column/rovv in which both of 
them have a nonsero entry. 

If i anu 3 are adjacent rows; and ±< j , 
then B ( i , j j = 1 . 

If k and 1 are adjacent columns; and k<( 1, 
then B(1 ,k) = 1 . 

i.e., the upper triangular portion of matrix B gives 
row-no jacencies , while the lower triangular portion gives 
column an jacencies. 

The Banclwiuth of a band-diagonal matrix is defined 

to be (height of the band + width of the band - 1 ). 

Vr 


gW - 


In the algorithm, we first permute rows to get a 
new row-sequence which tends to decrease the Bandwidth. Then 
we consider 1 the transposed matrix and try to find the new 
column sequence of the original matrix, after determining 
which column appears first* The desired column must satisfy 
the criterion Cl. 
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Criterion Cl ; 

We have with us, the new row-sequence. Let the first 
row in this sequence be R. Let P be uhe list of columns 
having nonzero entries in R, Let Q = P. 

If a row has nonzero entry in either of the columns 
given in P, ell other columns o.f that row which have nonzero 
entry in them are included in list Q. All rows arc considered 
and list Q is made. 

Now, consider the row-sequence in reverse order, 
and cancel out columns from Q if these have nonzero entry in 
that row. If at the end of consideration of one row , we have 
only one column m list Q, that is the desireu first column. 
If we have empty list Q left, then we consider row-sequence 
from tho beginning and the first nonzero entry in any of the 
columns of original Q determines the first column in new 
column sequence. 

If, after getting the new permuted matrix, Bandwidth 
is found to be larger than the maximum number of nonzero 
entries in either a row or a column, then we remove row's/ 
columns . 

We begin by removing one row first, and then one 
column, if necessary; and so on alternatively. When a column 
is to be removed, we consider the transpose matrix and remove 
a row, which is the same as removing a column of the original 


matrix . 
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If a row (R) is to bo removed to reduce the BW 
of the matrix, it should satisfy criterion C2. 

Criterion 02 ; 

The element of R at the right or left end of the 
band should be nonzero. All elements diagonally above the 
element of R at the lelt end of the band should be zeros. 
So also, should the elements, diagonally below the- element 
of R at the right end of the band, 

MAIN ROUTINE 

Step 1 j Given the random matrix A. 

Step 2 i Calculate row and column sums of A. MAXM = maximum 
of these row and column sums, IFF = 0. 

Step 3** Determine row and column adjacencies, and fill 
the entries of matrix B, properly. 

Step 4: Call Routine P. 

Step 5: The new- row sequence is obtained. Choose a column 
IR according to criterion Cl » 

Step 6: Transpose the matrices m and B. 

A<- A r 
B<e-B f . 

GO TO Routine P. 

Step 7: The new column sequence is also obtained. Permute 
rows and columns of original A as given by two new 

f 

sequences. PRINT the reduced Bandwidth matrix. 

Step 8; BW = Bandwidth of the new matrix. If BW^MAXM, GO 


TO Step 14 
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Step 

Step 

Step 

Step 

Step 

Step 

.No te : 

Step 

Step 

Step 

Step 

'' ' ' " 4 H 

Step 


9: If IFF is even, GO TO Step 11. 

10: A = A f . 

11: Check rows of A to find out which row satisfies 
criterion 02, so that it can be removed. 

12: Retranspose the matrix A, if it was transposed in 
Sxep 10. 

13: Make the removed row or column of A as "ZEPO" row 
or column. Call the new matrix A. Increment IFF by 
1 and GO TO Step 2. 

14: PRINT the final A matrix. STOP. 

ROUTINE P 

Whenever columns are being considered for resequencing, 
we are operating on A = A' (original), so treat sorted 
row (SR) vector as sorted column (SO) vector m that case. 

1 : Matrix A is available . 

2: Choose a row which has least number of rows 

adjacent to it. If there is only one such row, i.e. 
no clashes occur, GO TO Step 4. 

3: Choose the first of the clashing rows whose row-sum 

is minimum, 

4: Gall the chosen row IR. 

If degree ( IR) ^ 0, i.o. at least one row is 
adjacent to it, GO TO Step 6. 

5: Put Ik at the first unoccupied position from the 

end, of SR vector . GO TO Step 2. 
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Step 6j 

Step 7: 
Step 8; 
Step 9: 
Step 10: 

Step 11: 

Step 12: 

Stop 13: 

Step 14: 

Step 15: 
Step 16; 

Step 17: 


If we are considering columns of original 
matrix, GO TO Step 9. 

Put IR in SR VECTOR. If all rows over, GO TO Suep 24. 
GO TO Step 10. 

First column is already known; Call it IR. 

Modify the incidence matrix B for row IR. If I is 
incident to IR, put B(I,IR) or B(IR,l) = 0 depending 
on KIR or I>IR. 

Find out neighbouring rows of IR and put them in 
i'JGHBR vector. 

If NGHBR(1)^0 ac least one neighbour exists, 

GO TO S.ep14. 

Call the node following IR in SR- vector as IR. GO 
TO Step 10. 

Of the neighbours , select that row whose remaining 
degree (given by modified B matrix) is minimum. Let 
IR be that row. If there are no clashes , GO TO 
Step 20. 

If we are considering columns , GO TO Step 18. 

Choose row having maximum intersection with highest 
labelled node so far. If there are no clashes 
involved, GO TO Step 20. 

Of the clashing rows, find EXCLUSIVE ORS with 
previous labelled rows in reverse order. Whenever a 
minimum EX, OR is found, call that row IR. If no such 
minimum is found, call the first of the clashing rows 
IR. GO TO Step 20. 
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Step 18: Of the clashing columns, find out that colunn 

whose row-code appears first in already sorted-row 
"SR" vector. Call it IR. If there arc- no clashes, 

GO TO Step 20. 

Step 19: Of the clashing columns , find out one which has 
maximum row-code distance in already sortoi-row 
"SR" vector. 

Step 20: Put IR in SR- Vector. If all rows over GO TO Step 24. 

Step 21: Corapre-ss HGHBR vector removing IR from it. If 

there arc no rows in RGHBR vector, GO TO Step 23. 

Step 22*. Modify adjacency matrix B in so much that adjacency 
of IR with other rows pres unted in NGHBR vector only 
is affected. GO TO Step 12. 

Step 23: let DEGREE (IR) =0. GO TO Step 13. 

Step 24: RETURN to Main Routine. 


Example : 


Say, 


we are given a set of linear equations . 

Xj + a^ 2 Xq = 0 

a 21 X 1 + a 22 x 5 = 0 

a_ 1 Xj + a J2 x g =0 

a 41 x 5 + a 42 x 8 = 0 

a 5 1 x 2 + a 52 x 3 ^ a 53 x ($ ^ ^ 

a 61 x 4 + a 62 x 9 = 0 

®71 x 7 + a 72 x 9 


0 
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a 81 X 1 + a 82 x 5 + a 83 x 8 = 0 

a 9 1 X 2 a g2 , = 0 

The coefficient matrix A is ... (9 x 9) 


\ 


A B 


C D E F G H I 


a 



0 


0 0 


0 0 


1 0 


1 


b 

c 

d 

e 


1 0 0 0 1 

0 0 10 0 
0 0 0 0 1 

-01100 


0 0 0 0 
10 0 0 
0 0 10 
10 0 0 


f 

8 

h 

i 


0 0 0 1 0 0 0 0 
0 0 0 0 0 0 1 0 


1 

1 


1 -o 0 0 
0 10 0 


10 0 10 
0 10 0 0 


* 


When converted, tj a Bluck Diagonal matrix by row and column 
permutations , a becomes 



H 

A 

E 

I 

0 

D 

B 

E 

0 


h 

1 

1 

1 

0 

0 

0 

0 

0 

0 


d 

1 

0 

1 

0 

0 

0 

0 

0 

o 1 


b 

0 

1 


0 

0 

0 

0 

0 

0 


f 

0 

0 

0 

1 

0 

1 1 

0 

0 

0 


g 

0 

0 

0 

1 

1 

0 

0 

0 

0 


a 

0 

0 

0 

0 

1 

1 

0 

0 

0 


i 

0 

0 

0 

0 

0 

0 

1 

1 

0 


c 

0 

0 

0 

0 

0 

0 

0 

1 

1 


e 

0 

0 

0 

0 

0 

0 

1 

1 

1 



Looking at the coefficient matrix now, one can easily 


solve the given equations. 



CHAPTER 5 

GENER-sTABLE CLASbES Of MATRICES 


We have used the shift register technique for generation 
of structured Matrices , for convenience^ s ake T we redraw 
the figure here. 



Fie ax r - i 

q t til 

The (i+l)—^ row of the matrix is obtained fr^m r — 

row using the structural properties of the matrix. The allowed 

0 1 1 U X1S? clX 1 G * * * 

i) Shift left; 

ii) Shift right; 

iii) Count; and 

iv) Compare . 

An attempt is made to characterize the classes of 
matrices that can be generated using' shift register tech- 
nique and these four operations. It is based solely oh the 






49 


structural aspect uf matrices, and seems to be a little 
vague m nature. The mathematical characterization of such 
classes may be possible and can be attempted m future. 

Initially, assume that only operation allowed is 
"shift right", end a "zero" is appended at the leftmost bit 
position. In that case, 

r i+i •"= °> r i o r > 

O -h 

where, r 1+ . = (i+1)— row 

, implies concatenation 

1 ) = i— row with a rightmost bit clipped. 

The class of matrices that can be generated in this 

only those 

consists j± «±£ matrices in which ther^ appear bands in upper 
triangular portion of the matrix, the element in lower 
triangular portion being zeros. 

Similarly, if the operation allowed is "shift left" 
only, and a "zero" is appended at the rightmost bit position, 
then. . . 

r i+i = nil ’ 0 - 

“th. 

where, r^(1^ ) - ^ — row with a leftmost bit clipped. 

'And, all matrices having bands running; from N-B to 
S-W; and elements in S-E portion of the matrix being zeros , 
fall in the class of generatablc matrices. 

If instead of "zero" being appended, we append a "1" 
at the rightmost or leftmost bit position in cases where 
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operation allowed is "shift left" only and "shift right" 
only respectively, then the two generatable classes of 
matrices are respectively,... 

i) all matrices having bands running from N-E to S-f; 
with elements m S-E portion of the matrix being nonzeros; and 

li) all matrices having bands running from N-W to S-E, 
with elements in lower triangular portion being nonzeros . 

Next we consider the case when the only operation 
allowed is "shift right", but now either 0 or 1 can be 
appended at the leftmost bit position; i.e., 


i+1 


= 0 


r i ( V 


+ 1 


r ( 1 ) 
in r' 


where + indicates OR. 

In this case , all matrices having bands running from 
N-W to S-E fall into the class of generatable matrices . 

Similarly, if we have only "shift left" operation with 
either 0 or 1 appended at the rightmost bio position, then 
all matrices having bands running iiom N-E to 3-W fall into 
the class of generatable matrices. 

Let us no w allow two operations , viz., "Shift left" 
and "shift right", with "zero" appending only. Then... 


r i+l 


0 




-f 


r iOi 


), o 
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e .g. , 

Shift right 
Shift left 


r 1 - 


1 \0 0 X 1 1 \ 0 \1 

\\ \ \ v 

~ r 0 >l\o o \ i y 0 

/ i / /o o/i y / o o 

0 0 <M 1 <(o 0 0 

— >- 0 o o\i 1\0 0 


-> 0 0 0 0 \ 1 1 \ 0 

— 0 0 O'] 1^00 


To calcalate distance between r 1 ana a’j_ + -|s the following 
algorithm can be used. 

Step 1: DIST = 2* of cluster of 1 »s in r-. 

Step 2: If shift right and 1 of r^ is nonzero , GO TO 
Step 5 • 

Step 3: If shift left and 1^ of r is nonzero, GO TO Step 5. 
Step 4s GO TO Step 6. 

Step 5; DIST = DIST-1 
Step 6: STOP. 


This class includes all those matrices which have 
bands of 1 1 s running from either N-W to S-E or D-E to S-W; 
but all bands are symmetrical in their run, i.e. all change 
direction of orientation at the same row. 

In the preYious case, if we have appending of 1 *s 
instead of O’s; then... 


r ±+i = u nep + r fi)’ 1 

To calculate the distance between i 
same algorithm can be used by replacing the 


th 


tmd 

_ CCKTRAL LIBRARY 

term ’'iten^ro 1 ' 

Acc. No. 


2969 1 
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by "zero". 

This class comprises of all those matrices which 
have bands of zeros running from either N-W to S-h or N-E 
to S-W; but all of them being symmetrical in their run, 

Now, if along with two operation, we allow appending 
by both zero and nonzero; then.... 


i+1 


= 0 


r ± ( 1 r ) 


+ 1 


r i ( V 


+ r i ( 1 X ) , 0 


4 * 


r f h 


), i 


e.g. r 1 = 

R. — 0 
L _ i 
I - 0 
L - 1 
R - 1 
R - 1 
L - 0 

R - 0 implies shift right, append a zero. 

To calculate distance between r_ and r . , .. , the 

-L 1+ I 

following algorithm can be used. 

Step Is If "0" appended, GO TO Step 2. 

DIST = 2* No. of cluster of zeros.' 

If operation is R and 1 of r^ is 0 , DIST = DIST-1 . 

If operation is I and of r^ is 0, DIST = DIST-1. 



stop. 
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Step 2; DIST = 2* No. of cluster of I’s. 

If operation is R and 1 of r. is 1, DIST = DIST-1 . 

r l 

If operation is L and ^ of r ± is 1, DIST = DIST-1 . 

STOP. 

The class comprises of all those matrices which have 
bands of zeros and ones running symmetrically (i.e. changing 
direction of orientation at the some instant) down the matrix. 

let us now divert our attention to two other operations , 
count and compare. 

It seems possible to generate any matrix having some 
type of "blocks" in it using these two operations. All these 
matrices form the class of gcneratable matrices using these 
operations. The most general block matrix is shown in figure 2. 


? M Six 



PI PURE 2 


Knowing all s J s, ± T s and j's, the matrix can be 
generated by keeping truck of number of rows generated and 
changing the shift register contents when one row of blocks 


is over. 
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Mow, say we allow three operations, viz., count, 

compare and shift right with either "0" or "1" appending. 

s 

We shall also put one restriction that (i+1) — row is gene- 

■ +V. 

rated by using shift right operation on i — row. 

It is easy to see that this class comprises of multi- 
ple-band matrices running irom N-W to £-E as shown in Pig. 3. 



G-iven the first row and all s J s and f r s, the matrix 

"tj3n 

can be generated by shifting right; the i — row and appending 

s t 

a "0" or "1" as the case may be to get (i+l) — row. 

Similarly, ii we allow shift left with either "0" 
or "1” appending along with count and compare, the class of 
generatable matrices includes all multiple-band matrices 
with bands running from F-E to S-W. 

If we relax the restriction that (i+1) — row can be 
generated from i— row by either shift left or shift right 
operation vnly, it becomes almost impossible to characterize 
the class, based en some structural property. Hence, the 
case is not considered here. 


/ 




CHAPTER 6 


REPRESENTATION OE MATRICES 


In the previous chapter, we described the classes 
of matrices that can be generated using shift register 
technique. The characterization of these classes was based 
purely on the structural aspects of the matrix. 

We now present some representations of matrices. 

The research on these representations was undertaken for 
the following reasons, 

i) To characterize the classes of matrices, which can 
be generated using shift register technique, mathematically, 
so that it becomes simple to determine, to which class does 
a matrix belong. 

ii) Given a 'matrix , to determine whether it is a combi- 
nation of two or more of the structured matrices (i.e., 
obtainable from two structured matrices by AND, OR or 
EXCLUSIVE OR element-wise operations on them) # 

iii) To represent a given matrix completely and uniquely 
by some method. Storing the matrix in its new form may reduce 
the overall storage requirement of the matrix. 

Literature survey revealed the following three types 
of representations, viz., 

i) functional Representation; 

ii) Polynomial Representation and 

iii) Designation- number Representation !>]• 
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We investigated a fourth representation m this 
thesis, namely? 

iv) Walsh-f unction representation 

It should be mentioned that these represenations of 
matrices proved to be quite ineffective, in that, it was not 
possible to determine to which class a matrix belongs using 
either of them. However, they are possible representations 
and future- work may be carried out along these lines though 
it does not seem to be very promising. 

Functional Re present at ion s- 

Let us take an example. Say, a row is 

1 0 0 1 0 1 1 1 0 0 1 (n=t$ 

It can be represented byj~B,rJ ; where, 

r = column number of last nonzero entry in the row, 
and B is a function given by 

f k if entry in k column is 1 and there are 

4 * 

0‘s between k— column and r — column. 

J th 

£/ r ) _ <T -k if entry in k — column is 1 and there are 1 's 

\ ij ii. 

between r — and k — columns. 

' 0 otherwise. 

B is applied recursively to k to yield the given row. 

For the above example, r = 11 and B is defined as 


E 



B( 8 ) = -6 


57 


B( 1 1 ) = 8 
B( 4 ) = 1 


B( 8 ) = -6 B( 6 ) = 4 

all other B 's = 0 . 


Now, we can store |TB,r^j instead of the complete 
row. We can go a step further and after assuming xhe maximum 
number of nonzero entries in a row to be = p, redefine the 
function B as 


B( 1 ) = 8 B(2) = -6 B(3) = 4 

B( 4) = 1. 


That way, we require p storage cells for each row, 
and if the matrix is quite sparse, we reduce the total s corage 

p 

requirement quite a bit from N to p.F. However, core required 
when we store only the nonzero entries is p.N for the whole 
matrix, and in comparison we do not gain anything. 


Polynomial Representation: - 

The polynomial representation of a matrix can be 
explained best by an example 

1 1 0 0 0 0 0 

1 0 1 0 0 0 0 

0 0 1 1 0 0 0 

M = 

0 0 10 10 0 
0 0 0 C 0 1 0 

0 0 0 0 1 1 1 

0 0 0 0 0 0 1 
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View the column j as x J and write row polynomials 


of the 

matrix 

row 1 j 

1 

+ X 

2; 

1 

+ x 2 

3! 

x 2 

+ X' 

4s 

x 2 

4* X J 

5: 

x^ 


6: 

x 4 

1 

+ X' 

7s 

X 6 



Combine all row-polynomials to write the matrix- 
polynomial as follows. Whenever a row changes, change the 

> th 4* Vi 

sign ot the coefficient. Also view the 0 — power of x in j— 
row as (highest power of x in (j-1)— row + 1). 

For M, the matrix - polynomial is 


1 + x 


x 


4 7 8 11 

x^ + x ‘ + x - x 


13 ^ 19 

x + x 


- x 


24 


- x 


25 


26 , 33 

-x +x 


So, the sequence representing the matrix is 
1 1 -1 0 -1 0.-0 1 1 0 0 -1 0 -1 0 0 0 
0 0 1 0 0 0 0 -1 - 1-1 0 0 0 0 0 0 1 

The Sequence represents the matrix completely and 
uniquely. 

However, the sequence needs 34 words of computer 
storage which is much more than what would be required if we 
store only the nonzero entries. But if we are able to define 
row-polynomial in such a fashion that the final matrix- 
polynomial is "periodic", then we shall certainly reduce the 


i 
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storage requirements. 

It seems that this kind of representation would not 
help much m case of large matrices of order 50, or more. 


Designed; ion- number Representation ; - 

The concept of designation numbers is very well 
explained by Ledley £l<fj . 

Lr Tr 

let E be the given 2x2 matrix. Then, designation 
number of E (#^E) is the sequence of 0*& and 1 ’s obtained by 
placing all rows side-by-side. 

Define 2k variables A-j , Ag , . ..jA^, A^.^ , • * • ’^2k’ 


l k’ k+1 

lr* V* 

the first k for 2 rows and last k for 2 columns. 


# A i 


1 (2 1 ^ o's), (2 i- ^ 1 's)[ repeated to give a total 


sequence of length 2 


2k 


Then, express E in terms of corresponding elements 
in-^iA’s. The expression of E obtained in terms of A’s speci- 
fies the matrix completely and uniquely. 

10 0 1 
0 10 1 
10 0 0 
0 0 10 


e.g. 


E 
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E = 

1 

0 

0 

1 

0 

1 

0 

1 

1 

0 

0 

0 

0 

0 

1 

0 

A 1 = 

0 

1 

0 

1 

0 

1 

0 

1 

0 

1 

0 

1 

0 

1 

0 

1 

1! 

CM 

<3 

0 

0 

1 

1 

0 

0 

1 

1 

0 

0 

1 

1 

0 

0 

1 

1 

MM 

II 

0 

0 

0 

0 

1 

1 

1 

1 

0 

0 

0 

0 

1 

1 

1 

1 

A^ - 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

1 

1 


■6 — A., A_ A_ A, + A. A 0 A 7 A* + A- A,-, A™ A. + a^ A 0 A^ A. 

12 3 4 1 2 3 4 1^34 1234 

H~ A «| A^ A^ ^ A-j A^ ^3 A4 

or E = A^ A 2 A 3 + A^ A^ A^ + A^ A^ A g + A^ A g A^ A^ . 

i'or matrices in v/hich the final expression of E 
becomes simplified, it may be worth storing it instead of 
the full matrix. However, it is observed that for a random 
matrix, the expression does not simplify to a great extent 
and hence this representation is of no use as far as storage 
requirement reduction is concerned, 

Walsh-f unction Representation : — 

A system f(j,x) of real and almost everywhere 

nonvanishing functions is called orthogonal in the interval 

x -C. x x H if 
o— 1 

, X 1 

f f( d > x ) f(k,x) dx = Xj S 

■ ^ ^. k = 1 for j=k. 

= 0 for j 4 A. 
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If = 1; functions are said to be ortho normal Jf§ j. 
One of the sets of orthonormal functions is Walsh 
Functions. They are shown diagrammetically in Fig, 1. 

There is a close connection between f Sal f and f Sine 1 
functions as well as between 'Cal' and 'Cosine 1 functions, 

Wal (2i, 9) = Cal (i,0); Wal (2i-1,0) = Sal (i,0) 

i = 1 , 2, 3 , ... 

The product of two Walsh functions yields another Walsh 
function. 

Wal (h, 0) . Wal (k,6) = Wal (r^) 
where, r = h ® k in binary. 

e.g. Wal (6,0). Wal (12,0) = Wal (10,0). 

6-0110 
12= 1 100 
1010 =10 

Walsh functions form a commutative group with respect 
to multiplication. 

Plenty of such properties of Walsh and ^hlsh-Paley 
functions, and their applications are available in references 
|4, 10, 20, 58,44 ? 45 f 46,47] . 
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1’h.e first four Walsh elements are 


a .| - 1 1 1 1 

a 2 *= -1 -1 1 1 

= -1 1 1-1 

a 4 K 1 1 "1 



We stick to ternary (-1,0,1) number system and view 
the real line as 


» » 1 ! j j 1 -4 

- 3 - 2-101234 

0 1 -1 0 1-1 0 1 — 

= 1111 
# -1 -1 1 1 
- 1-1 00 

We also adopt a notation that if a row R = 1001; then 
R = + e .; i.e. e's specify position of 1 1 s in a row. 

( 1 ) gives us 


a 1 

= e 1 

4 

e 2 

4 

e 3 

4 

e 4 

a 2 

= -e 1 

Sm 

e 2 

4 

Qry 

5 

4 

^4 

CD 

a 3 

= -e 1 

4 

e 2 

4 

e 3 


e 4 

a + 

* e 1 

- 

e 2 

4 

6-z 

3 

- 

e 4 


-4 

e.g. a 1 - a g 


4 


Solving for e's gives,.. 
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e 1 = a l " 

°2 ~ a 1 " 

a~ — a. + 
3 1 

a 4 ~ a -i + 


a 2 

" 

8-^ 

D 

+ 

a 4 

a 2 

+ 

City 

5 

- 

a 4 

a 2 

+ 

a 3 

+ 

a 4 

a 2 

- 

a 3 

- 

a 4 


( 3 ) 


Now, any 4- column row can be expressed in terms of Walsh 


el or ont s . 


Q.g. 1 0 1 1 - e .j + + e^ 


= 3a 1 + a ^ 



a 


4 


a. 


a j + a 4 


for first 8 Walsh elements, the Table is shown in 
Figure 2. 



J! 1 

I 2 

P 

I 4 

P 

2L 

P 

- 4 8 


a 1 

1 

1 

1 

1 

1 

i 

i 

1 


a 2 - 


-1 

-1 

-1 

1 

i 

i 

1 


a^ -*■ 

-1 

-1 

1 

1 

1 

i 




a 4 "* 

1 

1 

-1 

“I 

1 

i 


' -i 


a 5 - 

1 

-1 

' -1 

1 

1 



i 

i 

a 6 

-1 

1 

1 

-1 

1 



i 


a ? - 

-1 

1 

-1 

1 

1 


i 



a 8 ^ 

1 

-1 

1 

-1 

1 


i 




FIGURE 2 
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e.g. e 6 - a 1 + a g + a 3 + - a 5 - a g - a ? - a Q 

a 4 = e 1 + °2 ~ e 3 “ e 4 + e 5 + e 6 " e 7 ~ e 8 
and so on. 

Now any 8 bit row can be expressed in terms of these 8 
Walsh elements. However, there is one snag. 

Say, w c are expressing N bit numbers by N Walsh- element s, 

N = 2 n , 

N 

If ^ leaves remainder 1 ; we get positive and true result, 

N 

as in the case when N = 4 or 16. But if j leaves remainder 2, 
we get negative result (*.* 2 on real line is -1 in our system). 

e. g. for N = 8 

— a^ -a g — + °a^ + a^ — a^ — a y + ag 

= -1 0000000 

N 

Hence, for representation of such rows whose — leaves 

3 

remainder 2, we negate the Walsh elements to correctly represent 
it, 

e.g. 1 0 1 1' © 0 1 0 = e 1 + + e^ + e ? 

= a 1 + a 2 + a 4 + a 6 - a Q 

True representation is -a. - a g - a^_ - a g + a g 

S im ilar tables as shown in FIGURE 2 can be developed for 

N - 16, 32 and so on. 

Let us now represent rows of a Band Diagonal matrix in 


% 


terms of Walsh elements* 


66 


1 1 1 0 0 0 0 
11110 0 0 
0 11110 0 
0 0 11110 
0 0 0 1 1 1 1 
0 0 0 0 1 1 1 
0 0 0 0 0 1 1 


Negate all results 
as ^ leaves 
remainder 2. 


R. 


R. 


IL 


(e . 4~ e o + e^ ) a™ — a. *t a*- — a. 


^1 t >=2 T °'3'' " "3 ” “/j. f “ °g + a Y “ a g 


( e 1 + e 2 + + e^) 

( e 2 + e 3 + e 4 + e 5^ 


a i + &2 


a. — a Q -f~ a„ ■{■ 4- Ui 


7 


R 


4 


a -| ~ a 3 


R 5 " 


a 2 + 


R 6 = 


“ a 2 


Rj 

% 


a^ + a^ + + a^ + a~ + ag 


a ^ 4“ a ^ " a~ — a^ . 


T/c. find that oven in terms of Walsh elements, there is 
no rovulority in expression of rows. And hence, storing waish- 
elemcnt equivalents instead of rows themselves does not lead 
to cor o-o.l location saving. 


Matrix polynomial was discussed in Polynomial represen- 
tation discussion. It consists of a series of I’s, -1 *s and 
0's, and hence that too can he expressed in terms of Walsh 
elements. As earlier representation, thou^i this is not very 
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helpful from storage reduction point of view, is a possible 
representation of matrix polynomial. 

In Chapter 2, we considered 6 structured matrices, not 
all of which directly appeal- in the application areas* How- 
ever, sometimes it is possible to express a seemingly random 
matrix in terms of two structured matrices, and core required 
to store two generatable matrices is less than that for a 
random one. 


g. 





1 

0 

0 

0 

0 

0 

0 










1 

1 

0 

0 

0 

0 

0 










0 

0 

1 

0 

0 

0 

0 







M 


















0 

0 

1 

1 

0 

0 

0 










0 

0 

0 

1 

1 

0 

0 










0 

0 

0 

1 

0 

0 

0 










0 

0 

0 

0 

0 . 

0 

0 





1 

1 

0 

0 

0 

0 

<f 



1 ' 


0 

0 

0 

°/ 

/.I 

1 

1 

0 

0 

0 

0 

0 



1 


N? 

0 

°y 

X 

1 

0 

0 

1 

1 

0 

0 

0 

A 


1 

\ 

1 

r 


A 

1 

1 

0 

0 

1 

1 

0 

0 

0 


0 s 


1 

1 

1 

V 

A 

0 

0 

'oi 

1 

1 

0 

0 



0 

cN 


1 

V 

X 

0 

0 

0 

0 

1 

1 

0 

0 



0 

0 



X 

0 

0 

0 

0 

0 

0 

0 

0 

0 . 



0 

0 

0 

0 

0 

0 

0 


Such representation is especially useful for large 
matrices. 

An algorithm was developed to find out the combination 
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of structured matrices, which would give the random matrix, 
but it is quite vague and hence is not mentioned here. 

In this chapter, we have described our unsuccesful 
attempts. The work along these lines may be carried out in 
future but does not look promising. 



CHAPTER 7 
CONCLUSIONS 


It is really surprising that though many authors 
have mentioned the occurrence of structured matrices in 
different fields, none of them has dealt with techniques of 
reducing computer storage using structural properties of 
such matrices. In this thesis, we have given techniques of 
generating, recognizing and storing structured Boolean 
matrices, with the intention to reduce the computer storage 
requirement of such matrices. 

Basically, we have considered six structured matrices, 
and have developed methods of generating these, using shift 
register technique. Apart from these we have presented 
algorithms for .... 

i) recognizing different structured matrices from their 
row or column permuted matrices, 

ii) recognizing three of these structured matrices from 
their row and column permuted matrices, and 

iii) reducing the Bandwidth of a given random matrix by 
possible removal of rows and columns. 

The generatable classes of matrices using shift regis- 
ter technique are characterized on the basis of their struc- 
tural properties. To characterize the classes mathematically, 
some representations of matrices e.g. Polynomial representa- 
tion, Walsh-f unction representation are given. However, these 
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proved to be our futile attempts and problem of character! 
aing the classes of matrices considered in this thesis, 

remains no solved* ‘ 



Q k» ww «a> 

O— * GENERATION OF BANDED! AGONAL mats V 

0 |l»>'» |S$» pm 

C— - N OT AT IONS ABX ' . . . 

•C— “ N = SIZE- OF THE MATRIX 
o-** ■ I =HOP I Z ANTAL DIMENSION OF TH- BAND 
C-—“ J=V OPTICAL 01 MANSION OF THE RAND 

jf" IK» HO* fc-.-S 

INTEGER • A( 30) ,R(3C) 

READ 101, N* I ,J 
PRINT 102, N, I,J 
PRINT 10? 


K = 1 

DO 5 10*?. , I 

5 A ( I D ) = 1 

IF ( I . £Q . X) GO TO 7 

11=1+2 

DO 6 I D»T 1 , N 

6 At ID)=0 

7 PRINT ICO, (At II ) , II = ‘ ,N) 

K = K+ 1 

IF ( K.GT.N) STOP 
DO 8 I D=3 , N 

8 B ( ID)=A( ID) 

IF (K.GT. J) GO TO IF . 


A ( 1 ) =1 

liito CONTINUE 

SHIFT RIGHT OPERATION ! DONE BELOW. . . 

DO 9 ID=2,N 
. 1 01 = 10-1 

t*|X9 ; . A ( I D ) = B ( I D 1 ) 

pliftici? GO TO 7 

iiiii A < I ) =0 
lligi GO TO 20 

fcilloS ? FORMAT (1 X,"*T HE SIZ'" OF THE MATR IX=* * 12/ 1 X , ®H OP I ZOI ' : A L 
■ '*' * 1 BAND DIMENSION **t I ?/ IX ,* VERT IC AL B AMD DI MANSION**;, ;IZJ 

i^JpX'FdRMAT t / 1 X , . *• THE GENERATED MATRIX .15. ».*) 

FOR MA T ( / 1 X* 3 0 12 ) 


'■ IHF COMPUTER- OUTPUT- FOR - THE J.£ST«DftTA IS... 

;XX : :\.'.. ; XX x ■ X; ■; X : ... X X.' * * ■>«- <&■ . w- *? 


■ TH6 SIZE OF'. THE VhaTRI'X : : . s l,0-\ 

HORIZONTAL BAND Of MlriSi ON- ” : 3, 
X VERTICAL BAND DIME MS I ON = «X g 


XX TH? G M M?T"^ IS.. . :lg;|x 







1 I 1; ! I 0. 0 0. 0 C 

1 l l ‘ 1 0 0 0 0 

0 ), 1 * ' } 0 0 G 

0 0 X 3 ) 1 .1 0 0 

0 0 0 1 1 t l l 1 ti 

0 0 0 0 I X 1 ’ ) 

0 0 0 0 0 1 l '■ 3 

0 0 0 0. 0 0 1 1 x 



GENERATION .OF BLOCK 


NOTATION 5 ': A?-' ... 

M = Sm OF- TH r MAT PI.: 
. M= NUMBER OF 0 LOCKS 
I, J = ARRAYS S p : C 1 F Y T 


OF TH r BLOCKS 


O'V: ' ;60 

m ' • 

IllliiiT 




TOT GFR ( : C) ,I( 0), J.C 2 ) 

PH AD lOlt NtMt.n IK), JIK),K=" ,M) 

PRINT S0:'»N,O 
PRINT 102 
L = T, 
t j = : 

L L =0 

ROW CRH ATI ON. DOMES HP TO : T AT 2 Mi’ NT NO. 8 

Do s io=:, tN . 

A ( 10=0 

IF ( L • £?Q • l ) GO TO 7 
LI = 1.-1 . 

DO fe 10=0 U 
LL=LL+I ( 10) 

JJ=LL+1 , 

N8=LL+nU 
DO 8 I D= J J » N8 
A l 10=1.. 

PRINTING THF BLOCK 
PRINT 100, (MTJ)fH-li ) 

L J =L J+ i 

IF (LJ.GT. J(U ) GO TO 60 
GO TO 9 
L = L + i 
LJ=1 : 

IFTL.GT.f4) GO TO TO 
.GO TO .1 

I L =0 

PRINTING REMAINING &ZERO& ROWS, IF. NECESSARY 

ID=i.H- . 

LL = LL + HIO 

IF ILL -GO. N) STOP 

LL=LL+1 

:|llif«L:.:GT.N) STOP 
liilsOliST D^IfN • 

. • ; 

lllp : :;I0 : 0 , ' ( ; M H ) , 11= 1 * O 

GO TO 1,22 2 

* E RHAO(/l>C>» *7011 GCWERaT'EP.'.HaTIITX TS.y.*) 

P M AT i n % » 3 01 2: ) ; o' 3 ggf §§ 

mmmmmsm i: 2 : 1 ■ ■ 2 . .... ' , 2- 2 '2 2 '■ ■ 2 : '2 -22'2': ■ -.22 ■ 
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. 2.0 - : 




« n » 


TH’" COM PUT 9 GUT put FOR T H TASTED*. TA 1.5 

THE Size 0 F T HP M ft TP I ~ =' 0 

T H f N l j MB' R nr Bt 1 r K S I M T H ' : ' m T R I 5 = * 

THF G- NfR'i. T r -0 Mft.TPlv. IS... 

1 0 0 0 0 0 0 0 0 0 . - 

I 0 0 0 C 0 0 0 0 0 

0 1 1 0 0 0 G 0 0 u 

011 0 COO 0 0 U 

o i i a o o. o o q u 

ooo] i : o o .o o .. 

0 0 0 0 0 0 1 a a. 0 : 

0 0 0 0 C 0 i ? 1 . 0 

0 0 0 0 0 0 1 1 a o • 

0 0 0 0 0 0 0 0 0 G 
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F 0 P M AT 

r I ; . 

CT: ) 



1 02 

F 

OR. HAT 

( /I v t 

► TH? ■' St Z"- : ! 

• i - 7 1 j r: c T y 

=*»I V 


■* 

TR ! ANv 

= 

♦ .T' ) 



1 03 

F 

HP MAT 

(/IK- 

* the G r o' :: ! 

!>/TFO HA TP IX 

IS.,,*) 

TOO 

F 

Q P m £ T 

i n x. 

3-CT? ) 




S'T 

ND 






£ ENTRY 

1. GOB o« o;c 


TH r CQMoyT R OUTPUT F Op TH TgST-DATA ' IS. . . 

j H f. si Z r - ■ 0 F T Hg M'A Tf? IX | S'. =10 
NO. /F TRIANGLES = ■ 

THE G' m ?-\ TFD MATS J X IS... 

0 o 0 0 0 0 

0 0 0 0 0 0 

€ 0 0 0 0 0 .: 

COO 0 0 0 

1 0 0.000 
1 1 0 0 0 0 • ' . 

1 1 1 0 0 0 . 

0 0 0 0 0 0 
0 0 0 0 0 0 
0 0 0 0 0 0 


1 0 0 0 
1 i o o 
0 0 1 0 
0 0 1 \ • 

0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 Q 0 
0 0 0 0 
0 0 0 0 



onoonoono 


— ■« ■' GE Ni; RATI ON OF l 1 i V R t A P P • *j G- B L Of. K ' 4 ?•. A T F 

Pi ,ti eca aw 

NOTATIONS AR-' .. . 

** — N-SIZE OF THE MATRIX 
«,«=«» m=no. of blocks 

— * I , J= ARRAYS GIVTNG SIZE '• OF SLOCKS 

K« ARRAY Gl y’H G FIRST COL UN M OF BLOCKS 

*«» 

INTEGER A { BO) » T ( 0), J( ; u) r K ( ” 0) 

READ 102* 0,M 

READ 101, ( It L) , JtL), K( L ) , 1. =•’ , M ) 

PRINT 1 0 , N» M 
PRINT 104 
L = 1 
LJ=1 

.1 CONTINUE . 

DO 5 I D=1. , N 
5; A ( ID ) =0 

NB=K( L) + 1 1 L )” 1 ' . • 

KA=K( L) 

DO 6 I D=KA , NB 

6 At I D ) = 1 

7 PRINT 100, (At II J,II = ) .,N) 

LJ=LJ+1 

IF (LJ.GT. J( U ) GO TO 60 
GO TO 7 
|p , L = L+1 
L J =1 

ISIS IFIL.GT.M) GO TO 70 
GO TO 1 

I iis* SZERO ROWS .GENERATED' IF - NECESSARY 

111® LL=0 

DO 10 I D=1 , M 
ll®# LL=LL+ J (ID) 
lliil TF(LL.EQ.N) STOP 
m m Ll=LL+I 
IpJlllS 111 IF (LL.GT.N) STOP 

111# -"'DO It:' IDal sN 

^K#:A'tID)=0 
■■lilt fitli§i#l 'iOO, (At 
t 

{$012}": ■ ' 

t / 1 X* » T. HE : s I ZF OF THE HAT P I.X =; 

lilllill- 

llllpiiiiv iiiiiiii©; iilll® |l • 

liiillilll 

ifiiiftr 



GO TO J. 
FORMAT 
FORMAT 
FORMAT 
. 1 BLOCKS 
FORMAT 


, I.|/ /'I 



th c 11 m p u t p output 


& 


T H T • 5 T DA TA 1 ■ » * * 


THE SIZE OF THE M TP /’ 1 S =' ; C 


NO. OF BLOCK r = 3 


THE 

G 

HI 

: *A 

T ! ‘ 

0 


, t ro | v i 

1 

\ 

l 

0 

a 

0 

0 

0 0 0 

. j 

% 


. 0 

c 

0 

0 

0 0 n 

1 

1 

% 

0 

0 

G 

0 

0 0 0 

0 

0 

1 

T 

1 


0 

0 0 0 

0 

0 

l 


T 

T 

0 

0 0 0 

0 

0 

0 

1. 

1 


r 

X 0 o 

0 

0 

0 

H 

jL 


i 

\ 0 .U 

0 

0 

0 

0 

0 

0 

0 

0 0 0 

0 

G 

0 

0 ' 

G 

G : 

0 

0 0 G 

0 

0 

0 

0 

0 

0 

0 

0 0 0 


GENERATION OF MULTIPLE 


n y : ' D | 


3 l :i c K. 5 & M A 


C« NOTATIONS. A r « « o 

O— N*SIZE OF THE M ATP T Y. 

C—« * NDV-NO • OF S- : TS OF BLOCKS 
C«— 5 ? F = ARRAYS GI-VlMG STATING 

C-— DIFFERENT BLO CK S . 

C-— N J =NO. OF BLOCKS IN A S T 


I NTEGER A ( 3 Q) * M (1 5 ) , S ( " ~ , 5 ) , F ( 
READ 101* N , I D V 
PRINT 105 * N,NDV 
PRINT 1,06. . 

DO 1 I=1,NDV 
READ 102, MCI ) 

KA=M ( I ) 

DO 2 J = i» KA. 

READ 103, SU 
CONTINUE. 

READ 104, (NJ (II ) »IJ a l»NOV) 

L = l 
LJ-1 

CONTINUE 
DO 5 I D=i, N 
M I D ) = 0 
KA=M ( L ) 

DO 6 K*l» K A 
N.BtK)»SiLf K» + F(L-,K.»-J 
IA=S(L,KJ 
IS-NB ( K) 

.00 7 ID=IA,IB 
: A(1D.)=1 

CONTINUE 

|| INT 100, ( A ( I . I ) t II = i ,-N ) 

I Bmm 

iFiLa.GT.NJa n go to 60 

I . |® TO 8i§ ' 8 

Ijl§§#+1. : 

*8I«lsi5T:i: : 

!lll||l!|f..«f^i GO TO 70. 


FINISHING. COLS. OF 


) t MB ( : ) , NJ( 10) 


iiiii 8 

88 T ; ;| 8 : 8 

ill! 


mam 

iiiips 


88 : 81g«i| : lT€' ; I.Zi«*.'''ROW.S . I;F :.NiCBSSA» Y' 

'■lilliliiiiii ■. 

L+NJUD) . 8 


p8s«S 

I 












X 02 

FORMAT (I?) 


103 

FORMAT (21?) 


1 04 

FORMAT (1012) 


105 

FORMAT (/lX,*THS : SJZ' 0- TH~ 

1 SETS OF BLOCKS =*»TTj 

s. s /. J r. T 

1 0 6 

FORMAT {/IX,* TH r; GFaTP tTFO M 

ATR I;: 

100 

FORMAT { /% X,3 01/ ) 

END 



SENTRY 

1003 

02 

0103 

0609 

03 

0103 

0405 

0709 


0207 

030204 


th5 computer output for th test-da ta 


THE SIZE OF THE MATRIX =10 
NO. OF SETS OF BLOCKS* 3 



THE GENERATED 

X X l 0 G 1 

I l 10.0 1 

i i; i o o i 

I I 1 1 l 1 
■ i a i i i i 

0 I 1 1 i x 
• 0 % 1 1 1 1 
0 X l l l 1 
0 1 11 . 1 l 
0 0 0 0 0 0 


M A JR I X . T S • * * 

I 1 ' 1 i 

111 1 
1 .1 ' l ' 

1 1 l 1 ■ 

X 1 J 
lio O' 

1 X o 0 
1 : 1 0 0 
1 1 0 0 
0 0 0 0 


o.— - 
o-~ 

O' 1 ™** 9 

Q «r« «® 

titmpv* 

C PPS(MR» 
V 

C‘ j *=••* 


GENFRATI DM OF V' S H ft P ™ D * ft TP I X 

NOTATIONS tip.- ... 

M = S I Z F OF TH C MATRTX 
n*MO. .OF COLUMNS In A ;T. DC OF \ 
Jl=MO. OF ROWS FORMING V- BLOCK 
A, B® ARRAYS * . WHICH WH-H C ONC 4 ^ 
DESIRED MATRIX- 


llli- 

jgjgfF| 


ife?-. 

mil; 

ISlIlili^ 






Wmmk 


Wm 


hhHh 


INTEGER -ft ( 1 r 'l ,8( 15 ) , C ( ; |,D( ’ P ) 

READ 101* • N,n. t J.l 
PRINT lOStNtUrJI 
PRINT 10" 

KK =0 
N J *N / 2 

IF (N«N1»R. EG. i ) GO TO ? 

K>N/2 . 

L*K 

GO TO 2 
K*N/2+l 
L = K“ 1 
CONTINUE 
M- 1 

DO 3 I = lt K 
A(I )=0 

IFn.LE.il) All) =1, 

CONTINUE 
DO 4 J = l® L 
Bl J ) =0 

IF ( J.GE.L-I1+I) 8 ( J ) = i 
. CONTINUE 

PRINT 100? { A { I ) » I = I » K ) * ( B ( J ) » J 
IF (M.6Q.N) STOP 
M is M + 1 

: v.. IF TM.GT.Jl ) RK=’ 

do Ws; i- 1 *k 

C(n«Ad) 

;Ii ^OO'tilC J = l» L 

D ( J)-Bt J) 

9aiI "f||il* ; K ; 

a ' 

F 'A ( I) -CtTI 3 

" =L* 1 ^ g. S 

Bl j") =D ( JO ) StS *»ft,: n ' ft" 

I 


Hi 


|» 






101 FORMAT (12) 

102 FORMAT (/1. 5',»TH'- SU K nf THF MATRIX =*,i ?./■ 1X» *N0. OF 

l COLS. IN HAL F‘ V =*■ , T : / T r *un. nf ft On? F 0° v! MG. V* BLOCK 

1 =* t 1 2 ) 

|03 FORMAT (/IX, » Th;: ^,Kvp >:-rto MAT R IV, IS...*) 

100 FORMAT (/l Cl? ) 
f?ND 

gentry 

10030? 


TH COM PUT; R OUTPUT FOR TH T ST*” DAT A IS... 

-T HE S I Z f D F T HE M ATR ! V 0 
NO. OF COLUMNS IN HAUF-V- = " ' 

NO. OF ROWS FORMING V-8L0CK •= - 

TH £ Gfc N f RAT" 0 M A Tp I X T S. . . 

ill o o o o i i 

1 1 l i 0. 0 /1 V I 

I 1 1, • \ . X ■ i 

1 i 1 t 11 l 'I 1 . . . . 

| 1 l V 1 • ) 1 l 1 

1, 1 1 T i 1 "I 1 '1 "* • . - ' 

I I' 1 1 1 V 1 l 1 ' 

0 1.11 1 1- 1 .1 'I 'O'.- 
0 0 1 1 1 11 l o 0 

0 0 0 1 l li 0 p 0 




Q IB* «SS# **■» 

C 0 N V ER 5 1 UN 'OF A G I V M M ?• f k ] J/'-jT ° PLOD iitT^JX FORM. 

Q* mas 

INTEGER • ( !*, 15) ,R ( T* ) ,cn. * ) f % F ; ( :■ *■ ) , SC( -».) ? ta GP f ]■ c ‘ ) t 
.1 B<15,15|, MSC (1") . 
q****#*NOTATI QnS » . * 

N sk 1 2, L UF- T H *M A i R T 
Qwhmhm tmvn A — I M PUT MATRIX 
jQwat-Bsn«8«e»«t!» pj =r ow sum v w» \ i , ip 
C ** *»f“*»*'C :s COLUMN $1 J M V P C T 0 R 
Qp»s.«»«MBRi2Wi f^j A XR t 1 M A XI M U.M OF i M • / ^ l ) M S 

c«p*-^«,maxc=makihum of th- column sums 

Ml NR ;=MINIMUM OF TH- ROW FllMS 
^ J NC- s MI N I M UH C F i H CO L L M j\i % y ( M S 
■—* S R ~ S OR T F D ROW VF CT'-ip 
•Q ns'* MM*** **» ^ q is SDR T E D C i )L UMN V ■? C ^ 0 R 

C#-**** 

READ -JLOltN 
PRINT IDE 
DO 1 , I ^ X f N 

. READ 103, ( All ,J>f J»T,N) 

PRINT 104, ( A( I, J), J=1 ,N1 



I II 


'0>V 

liisife 








WXSft K& «3S» WSS I 


CONTINUE 

COMPUTATI 0 N 0 F R 0 W SUM S A MO COL U M M S U M 5 

DO 3 I * 1 »N 
R (I )=0 
TAGitm-0 
DO 5 J»1,N 
P ( I )-R( I ) + A ( I » J ) 

DO 10 J*i,N 
C( J ) -o 

oo lo r«itN- 

| | ^ | | | | J' | 

•D6TSRMI NAT I ON OF MAXI MUN AND MINIMUM OF ROW AND COLUMN SUMS 

ill!!*®: 

M I NR = 20 
MAXC=0 

Mifi©*20 ■ ■ 

lliilii i=i>n 

I® I I jj|p{ M AX PU L T « R ( I ) ) 

SBpilf:ii.ST,.R{ I) ) 
lllliMAXC^Ll.CM) ) 

J||§|i!|||#I e CT : l > ? 

L i; ; cX*****' V CHFCK TO SE'f L^HETwER OTNDTTTON 'OF, BAND; DlAGONALTTT 

‘ ' ■■■ TiS : SATlSf I ! : : 0 *i; ■•■"■■■ 'L'L®L;|;: ; i ® 

■"'pHKIR,C,MA','.R,«^XC,M,IC5. L ■ L , V**®! 

\ rfj Trt i- ('- (1 

c; * " i . ; i f •»- * 4 , * ■■■ f * ) 1 j { n 

ISilllllllll^ 

WMMiiSkMBS0SiB9i 


:.R c l ; 
=R Cl) 
MAXC=CU 5 
MlNC-C I! I 


MAXF 
: Ml NR 



TC S = 0 

CALL R OWF I Ml M I NF , MA K.f , r , A f J ;« . j.\> 
IF (I.CS. a. i) GO TO - no 
CALL RWM A X R I M , M A K R , A T R , ' * Gk , Jf . 
IF ( ICS .FC. 1 ) GO TO *00 
IF (ICS. Rj. 2) GO TO ? 00 
C A LLR WDCRS ( A, R,T£GP , r C , ' » , MAXP * 

I F ( ICS * LQ..3 ) GO TO r 00 
)0 L = L+l 

I F { L • GT . if. ) G r I TO 1; 00 

IR D = MAXR" !i I : -JL 
MAXR=MAXC 
M I NR=MI NC 
DO 3,05 1*1 »N 
P. ( I )=C(I) 

T AGP ( I )-0 
SR (I ) = SCU ) 

DO 105 J*1,'N ■ 

35 B(I»JJ*AUtJ> 

DO 110 1*1, N 
DO no J = 1.,M .. 

L0 A t I » J ) =81 J , I ) 

GO TO A 6 
00 CONTINUE 
— REVERSE COLUMN ORDER IF NECLSSAR 
IF (MINC.E0. IRD+1 ) GO TO TIG 
DO 115 1=1, N 
15 MSC ( I ) =SC ( n 
DO 120 1 = 1, N 
J=N-I +1 

20 sem=MSC(j) 

10 CONTINUE 

^—PREPARATION for final printing 
I* do 125 1=1, N ; . 

DO 125 J=1 , N 
lilt. ,.jR =.SR ( I) 

IC*SC( J) 

|ii All, J) = BtIR,IC) 

ISIilsiiT«T 106;. 

Illilil.s : ii46 l*t,N 
lilKllili PtO*., ( A ( I , J ) , J=1 » NT 


PRINT 108 : 

' format , (id ;,; 

cnn M trr ; ( 1 1 « # f Rif ' G I VS ® '■ 

■11 ' 

1512:1 

iliilifiill 


MAllUX*'// 1 











IS ROUT IN!.' CHXK5 WH'.'H-R CC 
AGONAL IT Y IS SAT! 5FI r 5 qo mc 

S*= ♦ ICS = I, 


TcGER R(1S) f C ! ). 3 ) 

=maxr 

=MAXC 
LR =0 
LC=0 
LR1=0 
LCi = 0 
UNT1 =0 
IMT2=G 
?0 1= ♦ H 

(R ( I ) mi O.MA XR ) KOUNT* »K0UNT1-+} 

( C ( I ) . i.Q« Mft XC ) K0UMT?'K n UNT2+ * 

NT I HUE 

{ KOUNTi .NF. N^MAXP + I -0- . KOUNT? .NF . H-4 
(N-KOUNT! .! Q.Q) GO To ** 
(N-K0UNT1-SC. • ) IFLE = 

UNT*0 

=K2»1 

30 I =1 » N ' • 

(fill) .NE.K2 ) GO TO . 30 . . 

UNT-KOUNT+1 
IUNTI-KOUNTI +1 
INTINUE 

: ( I F L.R'1 • E Q «1 . AN 0 . KOU N T . Nf" . 1 1 GO. TO 5< 
'■( I FLR.EQ.G. AND .KOUNT. :Q.2) GO TO 25' 
(TFLR-EQ.1. ANO .KOUNT. -Q.l) GO TO 35 
: ( IFLR-fcQ.O* ANO.KOUNT. 0.1 > GO TO 31 

I TO 500 
; LR1=1 

) TO 25 ' ■ 

: (N-K0UNT2« C C.0> GO TO 45 
: ( N-KQUNT2.F C - 1 ) I F L C * 1 
HJNT=0 
5?=K3-1: 

I 40 I st ! j N 

tie (I }.*■!«. K3 ) GO TO. 4 0 
}UNT=K0UNT+1 . 

|li||^K0UN;?2 +l- - : T ;// 

JMf I'NUf • 

: ( IFLC I • e 0-1- ANO. KOUNT • N.S - 1 ) ' 'GO TO 5. 

l|||p: liofcit:Aii aK : OUNt i'ifel;T:' GO .TO tll 

|p || j:H||Sp ft NO . K,0U*4f . f'O .1;:| GO /.TCO 4,1 

?po:5bcv^:py;;::; 

= LG1 = I ' O.:::';;-: OIOOST’-- 

3 TO ; 3 5 : 

!|#!ttll:g:£ 3 so:,: oso: 

Wllllll 


THIS. R OlJ T I N' C HO S v S 0 W • WH H 3 r 3 r w* S U V 13 M 1 0 1 HU 1 ^ • 

INTEGER P ( }.E) ,SC ( 1'* J,T/. GP ( v ) 

DO 50 I = • » N 

IF (R (I ) . NF.MI NR.) on rn n • 

IP. = 1 

INCS'X— I NDFX+1 
SC ( INDEX) = I 
TAGR {.!)«■ 

RETURN 

50 CONTINUE 
END 

SUBROUTINE TnwFIN(Hl MR , 'ART , R f A., 7R » I ’’DR ,SC,TASR ,N, IC-S ) 

^wwra 

C— - THIS ROUTIN' FINDS" FIRST FT W ROWS n F MjB-TP. I X WHOSE' 

ROWrSUMS I NCR S AS FROM • ■ 1 1 NR TQ NftXR IN ST C PS OF )« IF SUCH 
C«>“* ROWS ARE NOT FOUND* ICS = I,' ELSE ICS-0 ’ 

^ wm tKK) 

INTEGER AUF, 15) ,R (IF ) , SCCT* 5) , T4GP <J3> 

55 IP “MI NR 

56 CONTINUE ' 

IFMP.EQ.MAXR ) -GO TO ?0 : 

DO 60 1= ,N 

IF(R(I).NE.TP+1) GO TO * 0 


K0UNT=0 



DO 65 J = I. ? N 

IF ( ( A ( IR * J ) • ; Q. O'. AND 
1 At I, JJ.EQ.O) ) KOUNT: 
CONTINUE 

IF (KOUNT.NE.l ) GO TO 
I P 3 1 P + 1 ■■ 

IR = I 

TN0ix=rNDrx+i 
Si t I NDEX J = I 
,.TA6Rm-l 
GO TO 56 
CONTINUE 
ICS = 1 
iRi TURN 

lsl|§ioo:l i me, wm&M t n » 


f’ JSfj 

. v- W 


MS: MHO-SI ';R0VI-SU'«»«.AXR.» AN® 
IF, .FOUND * »Cf= ! 0«. I.F NOT,,, t 


#i:|S RQUf‘1 NS. F INQS t N»«ftXR 
iiiliil ; FROM ; Pie V I OU S ROW' 

IF ALL ROWS OVER » ICS=2. 

INTEGER At 15, 1,5 ) *R t IF ) . 7<*GR li E) » SC 1 1 F ) 




^ . 


iiilSilllilHiiin'l 

IIS1H! 

■ ;■ 



DO fc 0 J= , N 

IF ( ( At IP : , J ) . : Q.O.Arin. &•( j t j j »pq. j .nr, . ( > ( IP, , j ) . AND. 

I A{ I, J ). ! Q. 0) ) KOUNT = Kf: .UNT+ 

80 CONTINUE 

I F ( KOUNT . N;‘ .? ) Gfl in 7 1 

K-K + i 
IR = I 

INDEX=INDFX+" 

SC ( INDEX) 55 1 
T AGR { I ) = 

IF ( INDEX. FQ. -I ) GO TO * 00 

GO TO 71 
75 CONTINUF 
ICS = i 
GO TO 85 
100 ICS=2 
85 RETURN 
FN C 

SUBROUTINE RWDCRS UGB ,'TAGR, SC, IP ,M AX* ,N ,1 NDFX) 

Q *R»f nan *** 

.£*»«••» THIS ROUT INF FINOS lOWtR PORTION. OF THE MATRIX WITH ROW?, 

C— — HAVING DECREASING ROW* SUMS IN 5T-PS OF J: AND DISTANCE B TwF 
C«—«“ CONSECUTIVE ROWS IF SuCH ROWS' NOT 'FOUND, I CS =1 f .'ELS- IC 


85 

86 


INTEGER A( I Ft 15 I ♦ R ( 1 5 » f 7 AGR ( 15) ,SC(i5) 

I P = M AXR*" J. 

CONTINUE 
DO 90 1=1 * N 

IF(R(I ) .NH. IP.OR.TAGR U J.Nf.O) GO TO 90 
KOUNT =0 
DO 95 J=1,N 

: IF s ( A( IRf J ) .E C. O.AND. A( I »J)*FC»l)«OR.(AlIR»J) . f C • 1 * AND . 
% A ( I , J ) . 0 . 0) ) KOUNT® KOUNT+1 
CONTINUE 

IF (KOUNT. NE.l ) GO TO 90 
TP=IP-I 
IR = I 

IN Cf:X = I NOiX+I 
iillllu SC I INDEX) =1 




ill 
|;)} ; : 

;; .95 

111181®!! 

iapi«® 

Ssfe,-®! 

Sill® 


i s li' ®'T;A6R ( I ) ' 

aili5 : '; ; !*FIINDfX.5,Q.N) GO TO 100 
: GO; TO Mb 

90 CONTINUE 

l®^;:iC : S®l , 

“lisiliitN 


000001,0 


THV COMPIJT P OUTPUT FOR TH T-- ST-OATA IS*. 

THF GIVEN MATRIX' 

1 0 0 7 C 0 1 

0 0 1 0 0 } 1 

1 1 0 0 1. 0 0 

0 0 1 0 0 0 
1. 0 0 • I - J 0 0 
0 0 .1- -1 0 0 1 

0 0 0 G G : : 0 

THE MAT P IK’ IN R'QUIR : D FOR M 

1 0 0 0 CO 0 

I 0 0 0 0 0 
1 i 1 0 0 0 0 
0 1 17 0 0 0 
0 0 1 1 l 0 0 

0 0 0 1 .1 1 0 

o ooo :i 


4 - 



Q*h«P»>Bv 

C'—~ CONVERSION OF a G1v~ 0 M :• TO I v I MI C:' A BLOCK Mi. TFT-; 

0“”“ ASSUMED THAT A COLUMN I : HOT COMMON TO ,vOR G THAO TWO BLOCKS 

TNT? GER A ( 1 f , | , B ( - ff * * •)»? (T 5 j , C 1 1 j . T , GR ( 1 5 | ,T & GC ( 1 5 > * 

l BLR < TO* ) » RLCt * 0 , ’• ) , MV< ? 0 ) , FLAG (TO) ,EB(! 0 >* SR P. f ) *$C ( V } 
INTEGER KTUO),F(lO,'tl) 

C*** NOTATIONS A R r A S F H t . L 0 W:, . 

C— - - A = INPUT MATS I X 

c *•*••“ b = dummy mat pi x 

R “ROW SUM VC STOP. 

C— “ C = COLUMN SUM V'CTOP 
C«— ■ T AGR=RQw TAGS • 

C*»-“ TAGCtCOLUMN. T AGS 

BLR=ROw FORMING A BLOCK 
Q «»* *»*» BLC = COLUMNS F ORM ING A P LOCK 
C- — . SB=S OR TFD BL 1 CK $ ' 

Ci™. ; SR = SORTED ROWS 
C»— - SC = SORTED- : COL UMNS 

jQ *83! essa ■» ' 


'• '' 

• ■ ■ 

1811 8 


rK 1 N t A U / 

L = C 
M = 0 
I CH = 0 
M1=0 

DO l 1=1* N 
R ( I ) — 0 
C ( I )=0 
TAGR ( I )*0 
TAGC(I 1=0 
read 1,03 » ( At- 1 * J ) »J=? *M> 
PRINT 104* (At I,J),J=1*N) 

calculate row and colum 

nn I si . n 



- » 



DO 53 I=: , L 
FL AG ( 1 )=(> 

KT ( I ) 

DO i>4 . J=: , L 
IF (I.£Q. J) GO TO 5-'-- 

IF IF ( I , J ) . EG. G) GO TO 
KT < I ) =KTC 
34 CONTINUE 
53 CONTINUE 


iT 


I B=Q 

DO 57 1= , L 

IF (FLAG! I ) » N~ »Q) GO TO 
IFtKTU >.GT.’ ) GO TO 5 00 
IF (Killl -GT. ) GO TO ' 7 
18 = 1 

CONTINUE 

IF (IB.EQ.O ) GO TO 500 
IB POINTS TO FIRST BLOCK 
INDEX 5 * I NOEX+V 
SB (INDEX) =18 
FLAG! I B ) =1 

IF ( INDEX.EQ.l- ) GO TO 3 00' 


iUCr. 5 


S B 


» i~ 


100 





CALL BLK0RD(IC»F,L,K7 } 

IF t IC.EQ.l) GO TO 55 

SB CO NT A INS BLOCKS I ? I OR D SR . 0 F P P I N T I N G 

INDEXR=0 

indexc=o 

IP =1 ' 

CALL 8LKARG( I NDFXR » INOEXC, I P* SB, L»8LC, SLR f TAGC, SC, 5 
CONTINUE 

PREPARATION FOR FINAL PRINTING. INCLUDE XZ SR 04 ROWS 
AND COLS, t IF REQUIRED. 






.DO .160 1=1 


IF (I.GT. INOFXR) GO TO .16! 
■:iR=SR.(i) 


mp£.. . .. „ 

: ©O' 170 J=1 5 *4 


I'F< J.GT. INDEXC ) GO TO 1*0 

IOSC(J) 


ma nit^f =b( ip, 1 c > 

:i!!««fiNUi ■. 
fera ltlii®Ns® 


:I F ( N. IQ.. INDFX.R.l .GO TO 175 

I' lDO '180 1 * 1 , m 


IX^l + INDEXR ■: ' ”■ 1 ;'vs : V;.Sv;^ : rivivi 





-..V- vii: -..vv vw v 
VllV ' : :v.5 : ; v;:;.v. : 

.■ .Vi 


u o w u 


DO 200 J = I, ,M 
■200 A' ( J j I Y ) =0 
195 CONTINUE • 

.190 CONTINUE 
PRINT :.o 
DO 205 1=1 t N 

205 PR INT I 0‘' *(A( I , J J » J 
STOP 

500 PRINT JOP 
■ STOP 

101 FORMAT 112) 

102 FORMAT ( ; .X t *T t* INPUT K '.TP I •;;«•■// ) 

103 FORMAT ( 5 T ? ) 

104 FORMAT C. 

106 FORMAT (//»!;<** THF- MATRIX' IN- P? CUTKr'O F OR M* » / / ) 

108 FORMAT ( 1 X HI 0 1 V c M MjTRJX CANNOT BF ' CON V CRT TO TO 
1 A BLOCK MATRIX*) 

END 

SUBROUTI NG 8LKF0R ( A, TAG- ,'-N, L, M t M/ = Rl.R f BtC, R , I CD 


THIS ROUTINE FORMS BLOCKS AND STORES RO^S AND C'US. 


F 0 R M I N 0 R L 0 C < 


0 ) GO TO 0 


INTEGER 4(1.5, 

N 1 =N‘“ l 

DO 10 

IF (TAGR 11 ) .N" 

L=L+1 

Ml =M 1 + 1 ■ 

BLR ( L j Mi ) = I 
■mGR(i)=i 
11=1+1 

lie IS J = II , N 
.M.l TAGR { J) .NS .0) GO TO 
J ) ) GO TO 

^iilrXt'O K=ltN 
IF (; ( I ,K)»A( J *K ) ,.N 
20 CONTINUE 

|It|l|St$ls» TAGR ( J ) = . 

DO K= ,M 
! WWtH-N 
liiil 


IN B LP ! I » . . i AND RICH ) RESPECT' 
5 ) ,TAGR ( l 5 ) , BLR (10*9) *BLC <10*5) »RM5) 


IlS 

n 

<W' 

i Si 


,0) GO TO 15 


1) SOTD'I 1 



Bill : ; 8 1 ■ i 


VFL Y - 



BLC,R,M) = K 
30 .CONTINUE 
M-G 
Ml =0 
GO TO 10 
35 ICH=0 
10 CONTI NUF 

RETURN 

fn c 


C 

C 

jQ «W 

c— - 


32 

SO 

45 

10 


SUBROUTINE 3LMmT(L,F ,Pt C) 

THIS ROUT IN- FINOS BLOCKS HAVING AT L~a:T 
BETWEEN THEM. IF BLOCKS ! ft. NO J T|T' : «5 r CT* 

I N TF.GPR FUG, 2.0 , BLC ! " 0 , " ) 

L1=L-1 
DO 40 I a1 » LT 
1 1 = 1+1 ' 

DO 45 J-I 1, L . 

DO 50 K=::,T 

IF ( BLC ! I * K ) » E C. 0 J GO TO 45 
DO 5 2 T 5 

IF I 8 L C ( J. t M).FC.O) GO TO TO 
I F ( B LC ( I * K ) . 0 F . BLC f J? M } ) GO TO 5 '? 

F { I , J)»l. 

F ( J ? I ) = 1 
GO TO 45 
CONTINUE 
CONTINUE 
CONTINUE 
CONTINUE 


A COL. IN COMMON 
FMtJ ) = F( J,I 


(J4*« 

• *-,VT4 - 

C *- , * w '' 1 ** 

!l! jjjjp|iT 


RETURN 

INC 

subroutine blkord(ICjF,l»kt» index, sb , flag, 18) 

; this routine arranges blocks in order in which they are 
Going to bf printed, ic-i implies a new block is to bf found 

AS PREVIOUS BLOCKS DO NOT INTERSECT WITH -.REMAINING BLOCKS. 

: 1C =2 IMPLIES BLOCKS ARE ARRANGED PROPERLY. 

liiliGfR; f( iG, iGitKiad) , ssiio l, flag 1 101 

r ■ 

Tlill® Ml it 

-6Q-C) ; GO TO -60 

' ■- 

I jlllT*^ ( J)|| | 

Illiu; i 





mok-^- ssoEi i GO : TO T OO ■ 


£10,1 j. Gi! TO O i 





CO TO 59 
60 CONTINUE 
TC = 

RETURN 
100 I C =2 

RETURN 

FND 

SU EROUT I ME BL KAR G( , I «C 


C"*" THIS R OU i I 
O- PROPER Of. OFF 


- A R R AN G c P.ijWS A*-'D f 
: FOR ■ F T ’-I/' L PR JOT] OR., 


f S P L * BLCtBLR » TAGC , SC., SR ) 
FROM SB- VECTOR IN 


INTEGER SR P 0 ) ,BLC < ' 0, 5 J , BLR ( . 0. , ' ) r T PGC ( i. F ) »SC (1-5 > ,SR (1: J 

99 K A = 5 8 ( I P ) 

IP1=IP+1 . 

IF ( IPl.EQ.L+l J GO. TO ? ; :: 

KB=SBIIPI) 

DO no 1 = 1,5 

IF(8LC(KA, D.cQ.O) 50 TO. .1.1 5 

DO 120 0=1,5 • • 

IF { 8 L C ( KA t I >. NF.BLC ! KB, J) ) GO TO X 20 
IZ-BLCl BA, I ) 

TAGC ( I Z ) =1 
GO TO 110 
12.0 CONTINUE 
iiO CONTINUE 
115 CONTINUE 

DO .125 1 = 1,5 

IFlBLCtKA, I J.FH = 0)+G0 Tpi. 130 
:iZ = BlC(KA, I) 

pp.. IF ( TAGC { I Z J .NE. 0 ) GO TO 1 25 . 

. INDEX C=lNDeXC +1 
. SC MNDfXC ) =IZ 

§1§1p continue ' 


I ■ '.CONTINUE 

mm f 3 t;H.s 

SM|ai PlF ( BLC ( KA, I ) * EQ» 01 GO . TO 140 
iMMft* 'iZ.«BLC{ KA» 1 1 . 

• IF ?TAGC( IZ ) .NE .1 ) GO TO 135. 

I NDf 1C. +1 

lHlll! SCXTNDf XC ) = I Z 

8BSp^ MlF?l»5' ' I »t 1.S ■ ■ 

1 11 * i ) life o i " sop ts :wq 

life- IP-IP + l 

poo P ; . ; : I F IIP . GT - L ).v G 0; ::T a- ,pPp‘ 

do p s: j - ■ » r i o " x ■ . : . fpv ..p : 

T m: tag r. t i t 1=2 ■ ■, ■ 







f-NC 

££NTPY 

07 

0000110 
101 0003 
ooooqoo 
0000001 
1010001 
0000001, 
0000110 


THt COMPuT. P OU T PUl FOP TH . T"ST-0*,TA- IS 
the Input matrix 


0 0 0 0 1 0 
1 0 1 0 0 0 1 • 

0 0 0 0 C 0 0 

0 0 0 0 0 0 i 

1 0 1 0 C 0 l 

0 0 0 0 C 0 1 

0 0 0 0 1 . 1 0 

THE MATRIX IN RFQUIP 0 FORM 


1 0 0 0 GO 0 
1 0 0 0 0 0 0 
1 l l 0 0 0 0 
l l 1, 0 0 0 0 
0 0 0 1 1 0 0 
0 0 0 1 .1. 0 0 
0 0 0 0 0 0 0 



o o o 


* ** « »>, Q t 0 kj y f : p 5 1 [ i hi n p £. 

Sf «K» IRiJ 

c#** NOTATIONS A* 5 "- A; 
£>-“ A« INPUT MATRIX 
C "***"" R = RQW SUM VfCTOS 
C-— • TAGRsROVi TAGS 
C— T AGC=COLUMN TAG! 
C SR s SORT SO ROW V* 

c— « p= Dummy vector 
O— B= CUMMY MATRIX 

INTEGER A { 1 5 5 3 5 
1 B (15 » 15) 

C~~- INITIALIZATION 

'Kiss «sa *577 

95 A4 lOljN 
7RINT 102 ■ 

DO 3 I=l; f N 
9<SI)=0 
T AGP ( I ) =0 
T1GCII )=0 
. R5 AD 103* (At!,J 
■ri TRINT 104, ( AH , 
9=0 

3AL3UL1TE ROw S. 
46 2 I =1 * N 
40 5 J = 1 , N 

lli: ' ft { I )=R ( I ) + A( 1 , J 

ili:- . MtRua.EQ.oj k 

IF(KR.EQ.O) GO 

| 4. J6NTIN4E 

PUT AZER04 ROWS 

KK=Q 

iitir i **i* n 

IF t ft ( I ) .NE.O) G 
ili KK=KK+I 
Tl sN-KK.ii ' 1 






INTfcRCHANGf C OLUMNS J ? ■ ^ T’in r ': 

DO 2.0 K= , M 
P ( K ) = A ( K 5 J ) 

A ( K » J ) = 4 ( K t I-IDEX) 

A ( K, INDEX) =P( K ) 

IF ( INDEX. FQ. M - K R ) 5 0 T n 1 0 0 
GO TO 25 • 

continue 

CONTINUE 
GO TO 500 
K = 2 

CONTINUE 

CHOOSE A ROW HAVING POW'-SU !i - K, A N D WH ! C. H I S N 0 T 
BET CONSIDER'-' D. 

00 30 i=; »N. 

IFtRf D.NE.K. GR.T4.GR l I ) ..NF. 0) GO TO 3 0 

K0UNT=0 

CALCULATE DISTANCE R r Tw‘ EN ROWS I AND V- 

DO 35 J*1 ,M 3 • 

IF { l A 1 1 1: J ) . EQ . 0. AND . A ( I P » J .} . FC. I ). OR . { A ( I « J) ■ F.Q* T . AND 
1 A (IR , J KFQ.O ) ) KOUNT = K: >UNT + t 
CONTINUE 

IF (KQUNT.NE.i > GO TO 30 
DO 40 J=l 1 N 

96 (Al I, J) .NEil.OR.TAGCt J) .NE. C) GO TO *0’ 

IF (AMR, J) .EQ.l > GO T Q 0 
TAGfU I)=i 
I N D f: X = I N D E X + J. 

TAGC( INDEX ) =1 
SR ( INDEX) =1 
IR = 1 
K = K+l 

INTERCHANGE COLUMNS L AND INDEX 
oci 45 L=1 » N 

P(L)=A( L, J ) 

MU J)=A!Lf INDEX ) 

A(L, INDIGO = PI L) 

IF ( I N D£X . F Q » N'-K R ) GO TO 100 
i GO TO 26 

IliiTlNUf 

lllIlliiUf- :; . 

liiitlie- ■■ -■■■ 

mMTINUE • 
iiiili :.imm ■" 

.33' 

vim 

' ; "-ACCORDING, i® , lisiiiiii : 


Sfi: imm 

sii'Mfi ; 

■ ' ? J f N . ' 3 , 3 3 3 ' 333' 33 ■ 3333/ : 


STOP 

500 PRINT 103 
STOP 

101 FORMAT (I?) 

102 FORMAT ( ! X * ‘■•T H r GIVTM M TR I v> : / / ) 

103 FORMAT (1511) 

104 FORMAT t IX* 15.12) 

105 FORMAT ( "> X * *T H? GI VtM- M'TRI: C A M NOT P-f . (.•DNVFR T r; D 
1 TO BAND T R I ft M Gli L A R F 0 R .v . # ) 

106 FORMAT (IX*// »* 1 H r M ATI- 1 X T N P '£l}7 c ''D FORM*, // j 

■ EN C 
£ENTRY 
07 

0100011 

101,0000 

000000 ! 

0100001 

1000000 

0000000 

0001000 


THF COMPUTER OUTPUT FOP THr TEST-DATA IS.. - 

THE GIVEN MATRIX 

0 1 0 0 0 1 1 

l 0 l 0 0 0 0 

0 0 0 0 0 0 1 

0.1 0 0 0 0 1' 

I 0 0 0 0 0 0 

0 0 0 0 0 0 0 

0 0 0 1 ; 0 0 0 


THE MATRIX IN REQUIRED FOR M 

JO 0 0 0 0 0 
1 1 0 0 0 0 0 
I- i ioooo 
0 0 0 1 0 0 0 
Olft 0 110 0 
0; 0 0 0 0 1 0 
1 0 0 0 0 0 0 



c~— PROGRAM TO CONVERT A G I V ” M .RAN 
C-— ' D I AGONAL FORM SUCH TH ft! a tl, J> 


BAND WIDTH OF THF CO.MV 

NOTATIONS USr.D ARF AS 
A- INPUT MATRIX 
N- SIZE OF THF MATRIX 
B- ROW COL . I NCI OF NC- 
FOR ROW AND L Owf R ON: 


T F ft R 


' F HR. TO 


a lo ws 


F ' |D COL. 


r-c D TP J 

mCIOFNC 


NGULAR PORT 


C— — F AND DUM- DUMMY MAT® IC S 
O— P. AMD C« ROW AMD COL. SUM VICTOR 


]Q im m» h» 

«K» 

^ tsssi K5a 
' C (ggi #m «Sfa 


DR AND DC- ROW AND. COL. DFGk T ur V r C TOP S 
SR AND SC- SORT! O- ROW AMD COL. VICTORS 
TAGR AND TAGS- '.TAG VICTORS FOP ROWS AMD 
ICL*» FLAG INDICATING WHETHER ROWS OR. COL 

BEING PERMUT'D. 


COLS, 


C-»-lFF- FLAG TO INDICATE' WHETHER ROW OR '.COL. IS BE I 
C-.— REMOVED FOR REDUCTION OF BAND*’ WIDTH. ^ v 


'j Q-fwa MS» nr*» 

j^«sa wa «*» 


M A X M- M A X I M U M MO. OF N 0 
I RW"® ROW OR COL. WHICH 
INTEGER A( 15* 1.5 > ? B( I r » I 


. 

O'** 


1 DR ( 15 ) , DC ( 1 5 ) , 5R( 15 ) * SC < 
INTEGER DUM! 05,13) 

READ 101, N 
PRINT 102. 

DO 1 1=1, N 

READ 103, ( AU »J ) , J = i,M ) 
PR INT 104, I A t- I, J ), J*'’ -»NJ 

INITIALIZATION PHASE 


BAND*’ width. 

-ZERO ENTRIES IN 
S TO 8C REMOVED. 
) ,FU*?** 5 ) , D UM ( 

1 1 *5 1 * TAGR (15 ) ,TA 


~P A R'i 


S R ( 1 5 ) , SC (15), TAGR < 


5 » T5 ) , 1 

GC ( 15 ) 


(IS) 


Qmummm 

IFF = 0 

IEE'dO'Z I-1,N 
llfiE; ' ;. DO Z J=1*N 

f. DU MM, J )“A( I? J) 
: :■% ' I€L=0 . .' 

00,5 t*l.N 

TAGC (ii — n 






0:1'.' E. 






1 10 J: 
; I ) =R c 
[ J)=C(. 






IF ( MAX M, IT.?; ( I ) ) MAX** = F (I J 
IF (MAXM.LT.Ct I ) ) M4TM »ct T ) 

12 CONTINUE 
IB CONTINUE 

CALL ROWING (A ,B,N»DP ) 

CALL COLINC(&,B,'i,nC) 

CALL NMBR W S i A , B , N , R , OR t r P # t ft, r,R , 7 Cl. , S C ) 
PRINT 10.7, (SR U ),! = ] ,N) 

DO 7 1*1, Ml- 
11 = 1 + 1. . 

DO 7 j = i: , N 
7 B(I,J)=B!J,?) 

CALL FRSC0L(5R,A ,N,TR ,7AGC) 

SC ( 1 )=IR 

CALL NMBCOLt N,F *0C» SC »T ACC, I CL , C , SR > 

PRINT 105 ’ 

DO 14 1=1, N 
DO 14 J =1 , N 


IC! ) 


14 DUNK I, J) = OUM(I, J) 

DO 15 1=1, N 
I R 1 = SR ( I ) 

DO 15 J=1 i N 
IC 1 = S C ( A ) 

OUMC-I, J)*DUMl ( IR 1, 

A(I,J) = FUR|, I.C1 ) 

DO 20 I =1 » N 

:‘tO PRINT 104, (At I,J ), J=! ,N) 

CALL RWCLRMtA ,MAXM,N, IR.W, IFF) 
■ IF(IRW.FQ.O) GO TO. 25 
9FF=IFF+1 

■96UFF«K*( IFF/2 5 .NF.,1 ) 60 TO 
DO BO J=i , N 
10 #(IRW,J)=0 

GO TO 3 , 
m 'CONTINUE 
DO 4Q 

40 A( J,IRW)*0 

Bill ■. 'GO' 'TO % 

llil; : continue 

■lliii 1=1, N 
Iiift5vva=/=l 

, J :) = DUM l I , J ) 


III 

III PRINT III " 
®1 11 .00R ; :S ; 0 : ::S'T":Si+i- 






STOP 

END 

SUBROUTINE PO*I.NC'( ft , R , rj , DP. ) 

C fiM pXL 

1 

C«— THIS ROUTlNf FINDS ROW TNCI D c ric~ S ft NO FILLS FNTRIES IN 

C""* UPPER TRIANGULAR P'1R T I0 ; ' ; OF B« 

Q Ml **» 

INTEGER .ft { l 5, 1,5 ) *B ( *• " , ' >, DR ( 1 3 ) 

N 1 -N**l 


46 5 I =T£51 

11=1+1 

DO 10 K , I 1 (N 

DO 1.5 J= : , M 

IF (At It J). R C. 0) GO TO 13 
IF (At I, JJ.NF. AtK,J )» GO TO 
DR ( I ) = DR( I H-3 
DR ( K } = DR{ K ) +1 
B ( I » K ) =1 
GO TO 10 
15 CONTINUE 
10 CONTINUE 

5 continue 

RETURN 

END 

2UBROUTIE5 IDLING ( A » B , N » DC ) 


. c— — THIS ROUTINE FINOS COL. INCIDENCES AMD FILLS ENTRIES IN 

lower triangular portion of b. 


hi ■ s 


■ INTEGER At 15» 15 > ».Bt 13 » 1 : > » DC ( 1$ )' 

3 3333; 

NI=N-1 

WmXm DO 3 J = l* Nt 
*S®331 Jl=j+1 

811: DO 10 L£J1*N 

iBSlli 3 00 315 I=lrM 

TfT At I » A ) . FQ. Oi GO TO 15 
lllli'3 . ^ I F ( A ( It All » L ) ) GO TO 1.5 

DC (L J “ OC( L ) +1 
lillli 3 3: DC ( J)=DCU)+I 
Simi 111! BCLf J 1*1:1 ■ 

GO. TO 10 


ill: CGMTINOt; 


30 CONTINUE 

■ return 

END 


■ r«: P ■■ ' 


i II l|§rll»l|:f illll I CLf,SSR ) 


' a ^ gor:1thm t 


p. NATION. U !; n IT '-TLiT 


NUMBERED 3 : 1 ; 


. : IM b n N ■ : A v* 

. A' ■ ■ .' ».** ^ /-* rthir t H F:‘ 0 


) , WL 


n s,, 

■ 3,- - : 3 iSltfSiiSS 


'''I-' . ' T ' '. ' '..3,. ■ 

33 33 ■ 3333 . 3;': 1 3T ' ' 1 3 ; - ,3 : 3 ■•:'■ ... : 3; ^ . :3 3 ..... -3:3. ■ ' 33.3 



O' — INhX- VICTOR CON T A Tn t MG NODES WHICH •'' Pr 
0»~ FINDING JNT'-R SECTIONS WITH HIGH". FT 

C- - 3 HR COL- NEIGHBOUR C?!. s. OF A COL. Ll.P c A. 

C ■ ARE BEING TEE T r D FOR fF at ro^COO- w ;; 

C--*» IN* SHOW- MP. OF M’* TGHB'HjP S OF A NODE B- 

c •***•■ ifl and i flF” fl ags 

C » 5D?9W- MF= OF P6W FSJ ?-H RPp-'RS; FIRST I. 

c ***••*• already found. 

INTI GER ■ US, 1 E ) * B U - 5 ! ' ) ,R ( X r ) , DR ( , - ) , ! 

I NGHBR ( V ) » T N rl'A ( If.'), ESP, H f ) , CCRC CL ( X F> ) , ■' 

KR =0 
! F LG=0 
I M C EX = 0 
I N =0 

DO 1 1=1, N 
j, NGHBR ( I ) ~0 

2 IFL=G 

CALL MNDGRW( T AGR ,OR,P , N , I FL, IP ) 


I F L = 1, IMPLIES ONLY OOF ROW OF M I N . DEGR; 

9F ( I FL = EH=A& GO TO 5 
.3 A L L MNR-WS M ( T AGP , R T N , I R ) 

CONTINUE 

IR CONTAINS R.OW NUMBER 


f..r. " 


9F LOR ( 91 ) =NF* 0) GO T T 10 
IZ»N-KR 
- SR I 1 Z ) — I R 
.T AGR ( IR )"1 
KR =KR + l 
GO TO 2 

lism 

GO TO 40 
CONTINUE 

IFTlCL.EG.il GO TO 6 

ifilNOiXeiNDIX+l ' 

iiiiiiiwotxi^Tt 

:a.: t agrt I R i**i 

Al-|ftL'|NBfl .fiO-N*-KRT HE TURN 

(CONTINUE ...'.A, . 

C A LL MODI F Y I T »B » JR ,10 , NGFSBRt 0% T AGR ) 

IFtNGHPRf? ) Q« 0 ) ■■'GO TO 40 

' df . 




CALL MNF v .O D ( M ,KDUNT» T H ! - ' , 1 MO = X,S?,P f ) 

GO TO , 0 

16 IR =NGHRR < 1 } 

! 5 INDEX- I NDcX + 

SR l INDEX I = ! R 
TAGR ( IR ) = ? 

I F U N C F. X . P Q . M • K 9 ) PF T U R N 
IFdN.FQ.lt GO TO TO 
DO 20 1= , TO 

I F < NGHBR { I ) . ! : » I R ) GO TO 20 
IP. 3 = 1 
GO TO 25 
20 CONTINUE . 

25 I R 1 = I R 1 + . 

■IF UR1.GT. IN) GO' TO. ’0 
DO 25 I = ?R).,IN 
I M = I -■ ] 

35 NGHBR ( I M ) = N GH BP \ I ) 

30 NGHBR ( IN ) = 0 
I N = I N**' I 

IF ( IN.FQ.0J GO. TO 3 9 

CALL MDFYP.W( INtOGHBR, 8,CR,IR) 

GO TO 12 
39 DR ( IR } =0 
M IR = SR 1 1 B 

I B=1 B+l 
GO TO 11 
m CONTINUE 

l«ir me 

immim PROGRAM COMES HERE only WHEN COLS. ARE. BEING PERMUTED 


m ma 


CALL FRSRWCf I N,NGH8Rf SS«> MIN* N, DR* A,NM8RW, JR» KOUNT , IK 
SLIT liORCOL, I MN ) . 

IF (KOUNT.FQ.l ) GO TO 50 

CALL MINRCOt I K* COR COL * IKN*NMBRW» 5SP * IR »N» A ) 

#8 CONTINUE 

plTggGG TO IS v. ■ 

iglfS* |ii ' 

Hill ; SuiROU T I ME MN DOR W { T AG R * €R * R » N * I F L , I R J 

RntiTTNP FINDS NODE' WHOSE .DEGREE I S MINIMUM. IF 

l^llllillfigOCCOR ? I FL =0. 

!!!!l7T;:Tit!la!2:’i *orti;5o *mis) 


MIN-20 




m-.Ti GO^TOgilli 

■TF (MIN- GT. DR Cl > 3 Ml««ORl I > 

S U ”n. : .N : . : __ i 

JFITAGR(i). 



g - 



IB l*ll« tMmm-y. 

1 ■« IFCTASfH H : > : -E.Q »T) GO TO I 
I F ( B ( IftSlf-T:O-.O)'--G0- : .TO l 


liii i itii 
<I) = DR( I )■ 
NT INI) 6 


■’ - ? • ;v 

i: : 1 B ■ - 

■B: : Bv : 


C*> • 

C" #w Ito* IS':' 

|Q 


IF =1 

CONTINUE 

IF (KOUNT.GT. ) p S' TURN 

! F L ~ 1 

R h TUR N 

MC 

SU PR 0 IJ T t N r, ¥ ' \ Fs W S M { T A G R , R , U , I P ) 

THIS ROUT I M r; FINDS OUT ROW WHOSE’ ROW* SUM IS MINIMUM 


INTEGER TAGS ( ).“> ,R O' ) 

M I NR *20 
DO 5 I=i r M 

I F IT AGR H } . R 0. . I ) GO TO/. 

IF (MI NR . GT „F { II) MlNF =F <7 > 

■ 5 CONTI NUF / ’ ■ . • B 

DO 10 I = " f N 

I F ( T AGR ( I ) . " : Q • I ) GO TO 0 . 

IF (R ( I ) .ME. MI -NR ) GO TO 0 

I R ~ I 

RETURN 

10 CONTINUE 

SUBROUTINE MODI F Y (N , B » I R* IN, NGHB R *DF ,T AGR ) 

C—*. THIS ROUTINE MODIFIES INCIDENCE- MATRIX OF ROW- NODE 
C»— THAT ROW WHICH IS -ALREADY NUMBERED . 


INTEGER B( 15, 15) ,NGHBR U 3 ) , DR (15 ) 
INTEGER T AGR I 15 ) 

DR (TR ) =0 

DO 5 1=1, .N 

IF (T AGR ( I ) . FQ .1 ) GO TO - 5 
IFU.6f.IR) GO TO 10 
IF ( B( I, IR) ..EQ .0 ) GO’ TO. 5 ' 

IN *1 N+l 
NGHBR( INI* I 
B ( I , I R ) *0 
DR ( I )*DRI I )-l 
CONTINUE 
CONTINUE 
lllilll: I;F;UR.1Q»N.) RETURN 

0 fe : F0«iBir4 •*» » n * m 


FOP 



O— THIS PC UT I N C H( 

c™ 9 ***** nf ighbou 1 s &\ 


■ ! j M ncf: } p 


■ pom ii khm G T HT 


Mr I b ri d U U b • T a L. p *•• t \ H Y h w g - ?••• 0 r 0 W • 

«P» KJWfcsa 

INTEGER M G HBP ( ]5 ) .OR l i f> 

I FLR =0 
M I N = 20 
DO 5 1=1, IN 
NR=NGHBR( I I 

IF (MIN.GT. DH ( NP ) ) Ml M -Of ( NR.) 

5 CONTINUE 

KOUNT =0 
DO 10 , IN 

NR=NGHBR(I ) 

I F ( K I N.-NF . DR{ NR ) ) GO TO 10 

KOUNT» KOUNT + 

IR=NR 

10 CONTINUE ' 

IF {KOUNT. GT. I ) RETURN 

IFIR=1 

RETURN 

END 

S U RROUT I MF I N T HU N( INDEX, SR * A, I N, NGH8R , M I N ».N » I R * K CUNT , I N LX 


in. CASE WHEN TWO 


HB'XJRS of .a row-nodf. have same 


O— - REMAINING DEGREE , THIS- ROUTINE. IS CALLED TO CALCULATE 
C— *? INTERSECTIONS OF THESE NEIGHBOURS WITH HIGHEST LABELLED 
£*»«•» NODE, AND SELECT IF POSSIBLE, ONE HAVING MAX. INTERSECTION 

Q*s« «©»««» 

INTEGER A( 15, 15 ) »SP (1 5 ) , NGHBP ( 15 ) , I 1 NTT 15 ) , NBR ( I *) , I NEXT 15 

INTEGER DR (15 ) 

IHLN=SR (INDEX ) 

INR = 0 

DO 1 1=1, N 
N BR ( I) *0 
INEX ( I ) =0 
1 iiNtm=o 

SBiiTi .DO 5 in* IN 

NR=NGHBR(I-J 

IFIMIN.NE. DR ( NP. ) ) GO TQ 5 

i a ■ 

iliil I NR ) =MR 
DO 10 |»,;N 

llllllllii iSiifil , J ) • EG- o ) GO T 0^ 1 0 

) . NE . At NR * j ) ) GO 'TO : m .. ' 

miii 

NTT NIL . ' 




KOUNT =0 

DO 20 !*•', im 

IF 1 1 INK I ) , M ;! . ma: 

kount=kgumt+ 

IR =NBR ( I ) 

iNfrX ( KOUMT ) = -i BR { 3 ) 

20 CONTINUE 
RETURN 
END 


T 0 ? 0 



SUBROUTINE 

K'XDRIM ,K TU 

;!T»lN tt X,TNDE> 

. ^ o R -f R 1 R f 

B3» Mft 

THIS ROUTT N 

c FINDS F v C L 

USIV OPS BE" 

ROWS 

«sa*S»3 

CLASHES AND 

ROWS ALRKaD 

Y NHJHPr R'-'D I' 

| p c yep j ;> 

£««•»«» 

WHENEVER A 

p j £ E ' J p D f- R 

OL VFD.» ROTUP- 

IS WITH MO 

Q CSD m* tt&i 

INTEGER ' INF 

X( IF ) ,EXOP ( 1 

)» SR UK, RE 

5 ) 


INTEGER Alif-t 15.) 

DO 1 1*1* KOUNT 

1 exoRm=o 
index 1= Index 

2 IQ=SR ( INDEX* ) 

DO 5 1=1, KOUMT 

IP = 1 NEK ( I ) 

K0UNT1=0 
DO 10 J=1 , N 

IF { A( I P f J ) »F0 .0 ) GO T 0 0 

IF ( A-( IP* J)*N£ *A ( IO»J ) ) GO TO 10 
KOUNT|=KOUNTi +1 
10 CONTINUE 

5 eXORU) =R ( IP) +R< IQ )-2*K0UNTl 

MN =20 

DO 15 1=1, KOUNT 

IF ( MN*GT. EX0R1I )) MN=EXOR ( I ) 

§im CONTINUE 

KOUNT 2=0 
DO 20 1=1, KOUNT 

Plfiftel IF TEXOR ( I ) .N“ .W? ) GO TO 20 

IR = J NFX ( I ) 

I K0UNT2=K0UNT2+l 
3&MTINUE 

I ) return 

iS5" : T5 v ';IN«X'l=.lNOEXl*K 

' llliliiiXK-W.. 0) GO TO 2: ; ■ ./K.TiBK 

li«l!ll« ■' 

TURN 

tiitlillf a? WTriR . ) 






DO 5 l = X, IN. 

NR=NGHBR ( I ) 

i F ( NG'HBR ( I ) . 5 T . I k ) G n T ? "t 1 0 
I F IB ( NR * I R J . C - 0 > GO TO 0 

DR f NR J =D^ ( ■ JR ) * I 
B l NR* IR ) - 0 
GO TO 5 
1.0 CONTINUE 

IF f B ( IR*NR ! . ’Q.O) G n iO r 
DR (NR)=DRC NR)*'- > 

B ( IR ? NR ) =0 
5 CONTINUE 


RETURN 

FNC 

SUBROUTINE FT SCOL ( SR *A»a*IP » T AGC I 

C— — THIS ROUTIN'” FINDS THAT COL. OF 1 WHICH WILL APPEAR 
FIR. ST IN FI NA L -M A If I •’ • 

Q «© W» 

INTEGER SR (1.5 ) »A(l5tT 5) » COLS ( 15 ) *T-AGC ( J i:; ) 


I A=0 

DO 1 1=1* N 


T COLS ( I ) =0 
IS : R=SR(M 
DO 5 v) = ItN 

IF (At I SR* J ) .1 Q. 0 ) GO TO . 5 



10 CONTINUE. 
:■ IA = I B 

: : 


'■‘‘ilT-”.?: 


Mr-- • •?**•>•••• •. wvmm 

laMi-il - 11 . ■■.I'll; 1; ■ 11 ,..l,i ill a v ill 

. . , ' • i- ' '•I-' ■ . ; : ; , , • ! ■' 1 . . ' ' 

«8aS»aO:5iii®':':E«' 

mmmrnm 

Iffllilfllaiasi^ 


ipipmi 


i*>=i A+i ■ • 

' a COLS ( I A)=<J 

'5a,alTA.GC.tu)-l :; 
i» sllla CONTINUE 
T : n : M=iA8. 

oiill 10 ill, n 

IF(I.Eti-lSR) 

1 -a:: ; leo 'Jf ill t m 

DO 5:0 K= ,N ' 

IFlA(ItK). 0. II GO in 

; «i r o 20 

I B=1 B+i 
COLStlSIi^ 

TAlC-TKllll., 

■ - 1 iiliitflltaiaii: ^1:: 

„.-j Z" • ' ! 1 . u .‘V 


TAGC ( J } =0 
30 CONTINUE 

IFlIA.EQ.l) GO TO '00 - 

IF(IA.FQ.O) GC TO 5 r . 

25 CONTINUE 

£.»<■>•» PROGRAM OEVTF C OmE S -HEP FR ON ' DOwNWAf DS. 

35 CONTINUE 
IQ =0 

ISR=SR ( 1 ) . 

DO 40 J= t N 

IF(A(I$R, J) » "C.Oj GO TO *0 


IQ=IQ+1 
COLS ( IQ ) = J 
TAGCUX 1 
40 CONTINUE 

IF (IQ.FC.l) GO TO 100 
I A = IQ 

DO 45 I ® 2 » N 
IST=SR ( I ) 

DO 50 J=? t I Q 
TC=COL$( J) 

I F ( A ( I ST? I C ) . NF . 0 ) GO TO 50 

T AGC (10 = 0 
I A* I A- 1 



f ^ F £ J f j| ) 

ALL. WM3R W$t A »B»N,C,nC, SCtTAGC, 1CL , SP ) 

Ml ioi /tsc m |j„ iv 

CjP.MAT i/lXt'- THE NEW COLUMN 5<QUcNCtr 


#, /IX, 1512) 


liiis 


MS 


isiiillll 


IFMA.6Q.1) GO TO 100 
CONTINUE 
CONTINUE 

PROGRAM NEVER COMES .-H$RF. FROM 45 DQWMWARD5. 

100 CONTINUE 

DO 55 I =.1 » N 

IF (TAGC (I ) . FQ.O) GO TO 55 
IR = I 
RETURN 
CONTINUE 

SUBROUTINE NMBCOL ( A, R*N,F»DG, SC, T AGC ,ICL,C, SR ) 

■ ■ 

iSKIil THIS ROUTINE TRANSPOSES A AND THEN CALLS NMBRWS TO 
^Ipl|#5 : ||NUM8ER ROWS (WHICH ARE COLS. OF ORIGINAL MATRIX) 

(l|ilS8ffii;5 : lW:|f Gf R M15, 1:5) , B (13 f 15 ) , F ( 1 5, IS) , S R ( 15 ) » SC( 15 ) , 
■illllllllii r.» DC ( 1 5 ) , T AG C (''1:5 I 



<>-«• FIRST IK SOFT SO ROW V CTHP, 

jp HP «S liS'Sl 

I NT* GfP. ijGHBP { 1. f } , S5 C > { i ' ) t DP ( IF ) *A (I 5,1 ' ) , 

1 NMBRWl r ) fCOPCOL ( * - ) 

I K =0 

00 1 1 = 1, N 
COR COL £ J > = 0 
1 NMBRWtl ) =0 
DO 5 1 = 1, IN 
NP = JGHB.n ( I i 

I F ( MI N . N? . . DPI NR) ) GO TO F 
DO TO J= , 0 
9SR = SSR.{ J ) 

IF(A(NR,ISR).N c .l) GO TO 10 
I K =1 K+3 
NMBRWt XK)=J 
COR COLI I K ) = NR 
10 ' CONTINUE 
5 CONTINUE 
IMM=20 

. DO 15. 1=1, IK 

IF I IMN.GT.NMBRWI I ) ) I MN=NMBR Wl I ) 

'15 CONTINUE 


KQUNI-0 
DO 20 1=1 , IK 

■ IF-iNHBRwi I ).NE. I MM) GO TO 20 
IR = COR COL ( I ) 

■K-OUNT=KOUNT+X 

MQ' continue 

RETURN 

END 

SU BROUTINE MI NRCDC IK,CORC0L, I MN, NHBR W, S-SR, IR , N* A ) 


pss# 

.£*' Es- m 

c*»* 


THIS. ROUTINE CALCULATES MAX. ROW-CODE DI STANCE BETWEEN 
s CttLS* WHICH HAVE CLASHED ON. FIRST ROW-CODE APPEARANCE- COUNT » 


. iHTiClR. ;CORCOLI 15 ) ,MMBRWU5) ,'SSR (15 ) ,A { 15, 

MM : 1 DIFF( ill , PIN UM ( 1 5 ) 

I IT=0 











M I N = 2 U 

DO IS I», f IT 

IF ( KIM.GT . DIF FI I ) ) HT N = CT FF l I } 

1,5 CONTINUE 

DO 10 1=0, IT 

IF (DIFFU ) .o: .MIN) GO TP "'0 
IF=RFNUM(I! . 

RETURN 

20 CONTINUE 

F N C 

SUBROUTINE RwCLRM { A, MAXfN N, IRw, IFF) 

Q KM) IBS* WiSf 

C~— “ THIS ROUTINE CALCULATES BAND- WIDTH OF A AW D CALLS 

c *»■»»» which rfnovfs rows and cols, of & to reduce: its ew 

Q TO> KTO 

INTEGER ‘ M l 5, 15 ) t D( 5.5 , 1 5 ).*BW ’ ... 

I R W = 0 
. NCLFR=0 
NRWFC=0 


C-— 

Q BiSt m BE® 

C* 


NCLFR“ no. of cols, in first row 
NR wFO NO. OF ROWS IM FIRST COL. 


DO S' 
DO 10 


N 

vKl* I 


,10 




T: 

I# 


1 v 

ffy- 


lilt 




IF ( At 1 , J ) * EG. 0 ) GO TO 10 

ifr.*i- j+i- 

IFUFR.GT. NR W FC ) NR WF'C * I FR 
GO TO 5 x 

CONTINUE 

continue 

' N1=N«1 
DO 1.5 1=1, Nl 
DO 20 J = I, N 
J J =N«* ( J"I ) 

I F ( A M , J J ) . o ) GO TO 00 

T | .25. 1 1 T 1 

. IF ULS.GtCnCLFR ) NCLFR-ILS 

E . .GO TO 15 
, CONTINUE 
fl I CONTINUE 
11 llltiiiFC^NCLFRrl' 

■mb bw ■. 

IF ( BW./6Q. MAKM ) RETURN 


IR J ~NR W.FC- 


i;T 


. TCT^NCLFf: 

I S*'#. T ' 


r it » i &mmmi t - .m. tpi . X r.l« T R w . 





V'E 


i OUT IN 



DO A 5 J= , N 
35 A ( I » J ) =D( J i T ) 

CALL RWRMVLC A *N t .IC1 , T F ’ . I R w ? IFF) 

DH A 0 1 = ,M 
DO 40 J= , N 
40 AU,J) = D!lfJ> 

101 FORMAT { /1X»* Trtl RAM D“ **, I DTK *♦,!?), 

R FT (JRN 
PNC 

SU 8R OUT I OF *?W RMVL ( & , 0 , I 1 , 1 CT , TR W» IFF ) 

C”"“ THIS ROUTT N- r REMOVES ■ ROWS /COLS . CF A TO REDUCE 
£*« w»a ITS BA ND ; - :: W I DT H » 

^ *Bffl e;s» 

I HUGER A! IS 15) 

I Q = 1 R I + 1 C 1 

DO 5 I =1? IRA 

I J = 1+1 Cl 1 

IF(IJ.GT.N) U=N 

IF ( At I* IJ) «N ::; .1 ) GO TO A 

IF M J.FQ.N) GO TO 20 


f 0 I N 


) GO TO 


25 

20 


liSIll 

Hi;r|A A- A 

AAi :AA';v AAA, 
4'AA : \A'A' ' ' 
sAA/'A A"; 1 
AAA'.:A:::: ; v',A' 

Gy All' ; _ ;i ’i‘ 3 -. ■ 


jIN=n-1j 
X 3L — I 
DO 25 K 
11 = 11 + 1 ' 

JX=I J+K 
IF LA ( II, Jl ) .Ft 

CONTINUE 
IRW = I 
GO TO 30 
CONTINUE. 

IF ( !R I.eO. M) GO TO 3 0 
IF ( 1R i+lCl . GE .N ) GO TO 
IF IA( IRliI) lNE.0) GO TO 
IRf.= IRl 

IC2 = IRI + 1'C!“I. 

IK -I 

DO 45 I -IQf Ml 
"I^ 2 *IR 2 +I 
li*«lC 2 +t- 
sliitt+tA:- :} 
illil IRlj IKKNS^I 
AiWlilRI' :A1 

i c i 8 i c i + 

i# j 8 1 c li 



OR. AIIR2+IC2) . NP 


% ) GO TO 


’ 50 




mil in it ii ii 

■A ; S'.:A ■ . ' . 



rONTINU! 

IP ?. =N*“ I C 

IFUQ.'G '.NJ ! P 7 ~ I R ** 

IC’=N + l-]Ri-” IC?. 

I F ( IQ. G~. .N ) 1 C2 = i ' 

I P 2- ft C . n F 3 0 W5 A L R r A D Y C flM $ ? 0 r Rr D . 
I C ?- N 0 . ' 0 F C 01 $ . TO B r LEFT MOW. 


I F ( I R 2 . EG . M J GO TO TO 
TD =0 

CO 60 I = IR 1 » I R2 

I D = I D+ 1 , 

IF U( I? ID) .N" .0) GO tq o 
60 CONTINUE 
I G = N*»- 1 R 2. 

DO TO I-* 1 » T G- 
IR 2 =IR 2 + : 

I C2= I C2+1 

IF(A(IR2»IC2I .NF.l.OR .AUR2,N} .NF.l ) GO TO 75 
I R W = I R 2 
GO i Q 30 
75 ' CONTINUE 

IF ( A- ( I R 2 » IC?) »Nf: . 0 ) GO 7 0 30 
70 CONTINUE 
30 CONTINUE 

IF < (IFF-( IFF/2)*2).EQ.l) GO TO 8C 
PRINT IOi, IRW 
RETURN 

80 PRINT TOE, IR.W 

1 01 FORMAT ( / / 1 X » *■ THE REMOVED ROW IS#? 12) 

|02 FORMAT THE REMOVED . COLUMN I.S®fI2) 

RETURN 
END ■ 

f ENTRY 

0100000 
001000 0 
Iltlill 


iiliil o 


mmmm: 




, . :7m 'COMPUTER: OUTPUT FOR' THE TEST 

|>|i ; | ' . | .7 : v; '3 

, ■: ■ . i ■■ , .'■ ■■ ■ 1 ■ ;■ . ■ 
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« £ L..n nSM asffiraSSJffi 
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